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These lectures, presented at the International School of Physics "Enrico Fermi," 
deal with two major themes. The first is the remarkable story of the solar neu- 
trino problem, which (along with the atmospheric neutrino anomaly) recently led 
to the discovery of massive neutrinos and neutrino oscillations, physics beyond the 
standard model. I will describe the physics of the standard solar model (SSM), 
the experimental program that was motivated by the discrepancies between SSM 
predictions and the initial observations of Raymond Davis, Jr., and his colleagues, 
and the recent results of SNO and SuperKamiokande. These first lectures end 
with a description of what we have learned about neutrino oscillations and the 
neutrino mass matrix, as well as the open questions (neutrino charge conjugation 
properties, the absolute mass scale, CP violation) that could ultimately impact 
our understanding of baryogenesis, the origin of large-scale structure, and other 
topics in cosmology and astrophysics. The second theme is the core-collapse super- 
nova mechanism and associated nucleosynthesis. This problem connects neutrino 
physics, which controls much of the nuclear physics of the star, with the long-term 
chemical evolution of our galaxy. In particular, the r-process, which produces 
about half of the heavy elements, remains poorly understood, despite important 
new constraints from studies of metal-poor halo stars. The possible role of new 
neutrino properties on both the explosion mechanism and nucleosynthesis is noted. 



1 Introduction 

These lectures were presented at the International School of Physics "Enrico 
Fermi," August 6-16, 2002, "From Nuclei and their Constituents to Stars." 
This written version includes a few new results from later in 2002, such as the 
recent announcement by the KamLAND collaboration. 

The main theme of these lectures is the interplay between neutrino prop- 
erties — their mass, mixing, and behavior under charge conjugation and CP 
— and astrophysical phenomena. The first topic is the solar neutrino problem, 
in which the discrepancy between the predictions of the standard solar model 
(SSM) and the results of the chlorine experiment ultimately led to the discov- 
ery of neutrino oscillations by the SuperKamiokande and SNO collaborations. 
The lectures include a discussion of the basic physics of the SSM, the exper- 
iments on solar and atmospheric neutrinos, the effects of matter on neutrino 
oscillations, and the current status of our efforts to determine the neutrino 
mass matrix. 

The second topic is one reminiscent of the solar neutrino problem in the 
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1970s, the core-collapse supernova mechanism. The failure thusfar to develop 
a robust theory of supernova explosions may indicate some basic inadequacy 
in our treatment of the nuclear astrophysics (such as our inability to model the 
hydrodynamics and transport realistically in three dimensions). But surprises 
could arise from new neutrino physics, as the environment (intense neutrino 
fluxes, high- velocity matter flow) is significantly different from any other we 
have probed. The supernova mechanism is important to basic astrophysics, 
as it controls much of the long-term chemical evolution of our galaxy: su- 
pernovae both synthesize and eject new elements. Futhermore, as a source of 
measureable fluxes of neutrinos of all flavors, supernovae offer experimentalists 
new opportunities to test neutrino properties. The lectures deal with super- 
nova neutrino physics and with the nucleosynthesis associated directly (the 
z/-proccss) and indirectly (the r-process) with neutrinos. 

The "live" audience for these lectures was advanced graduate students and 
postdocs in nuclear and particle physics: the material is covered at this level 
and at a depth appropriate to a survey. 

2 Solar Neutrinos 1 

The neutrino has been with us since Wolfgang Pauli's proposal, in Decem- 
ber, 1930, that the emission of an unobserved spin-1/2 neutral particle might 
explain the apparent lack of energy conservation in nuclear beta decay 

(A, Z) -> (A, Z - 1) + e+ + v. (1) 

Enrico Fermi was present at a number of Pauli's presentations and discussed 
the neutrino with him on these occasions. In 1934, following closely Chad- 
wick's discovery of the neutron, Fermi proposed a theory of beta decay based 
on Dirac's description of electromagnetic interactions, but with weak currents 
interacting at a point, rather than at long distance through the electromag- 
netic field. Beta decay was descibed as a proton decaying to a neutron, a 
phenomenon energetically possible because of nuclear binding energies, with 
the emission of a positron and neutrino. Apart from the absence of parity 
violation, which awaited discovery until 1957, Fermi's description is a correct 
low-energy approximation to our current standard model of weak interactions. 

The neutrino was connected early on to astrophysics. As Bcthe, Critch- 
field, and others unraveled the stellar processes for hydrogen burning (the pp 
chain and CNO cycles), stars were recognized to be copious sources of neutrinos 

4p -> 4 Hc + 2c+ + 2v c . (2) 

The detectability of neutrinos - of which Pauli had apologetically dispaired - 
was established by Cowan and Rcincs in 1956, and the existence of more than 
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one flavor of neutrino by Lederman, Schwarz, and Steinberger in 1962. Thus, 
with the development of the Glashow-Wcinbcrg-Salam electroweak standard 
model and its prediction of weak neutral currents, it was recognized that su- 
pernovae would emit prodigous numbers of neutrinos of all flavors: the hot 
protoneutron star at the core of the collapse contains a thermal sea of trapped 
neutrinos produced via neutal currents. Questions of astrophysics and of neu- 
trino properties are thus interwined, and much can be learned about one by 
improving our understanding of the other. 

More than three decades ago Ray Davis, Jr. and his collaborators 2 con- 
structed a 0.615 kiloton C2CI4 radiochemical solar neutrino detector in the 
Homestake Gold Mine, one mile beneath Lead, South Dakota. This experi- 
ment began a field, neutrino astrophysics, that in the past few years has pro- 
duced remarkable discoveries showing that our current standard electroweak 
model is incomplete. The new neutrino properties that have been established 
- neutrinos are massive and neutrino flavors strongly mix - provide our first 
glimpse of new interactions that likely reside at energies far beyond the limits 
of current accelerators. Furthermore, the discoveries have important impli- 
cations for particle dark matter, large-scale structure, and models where the 
baryon number asymmetry is connected to leptogenesis. 

The first few years of the Davis experiment showed that the number of 
neutrinos detected was considerably below the predictions of the SSM, that 
is, the standard theory of main sequence stellar evolution. Those results were 
refined over the 30-year lifetime of the experiment, which ultimately achieved 
a precision equivalent to about 3% of the SSM prediction. The CI experiment 
was followed by the SAGE 3 and GALLEX/GNO 4,5 gallium experiments, the 
Kamiokande 6 and SuperKamiokande 7 water Cerenkov detectors, and the SNO 
heavy water Cerenkov detector 8 . Furthermore the explanation for the discrep- 
ancy that was inexorably forced on us - massive neutrinos oscillating as they 
travel from the sun's core to the earth - proved to account for a similar anomaly 
in atmospheric neutrino measurements 9 . These new neutrino properties are 
now being confirmed with accelerator and reactor neutrino experiments, such 
as K2K and KamLAND. Ultimately a series of long- and short-baseline accel- 
erator/reactor experiments coupled with new astrophysical measurements will 
yield precise values for the mixing parameters and probe new phenoma, such 
as CP violation. 

My goal in the first half of these lectures is to summarize the solar neutrino 
problem - the standard solar model, the observations, and the accumulated 
evidence that led to neutrino oscillations as a resolution. The status of our 
knowledge of the neutrino mass matrix will be given, as well as the exciting 
set of outstanding questions - particle-antiparticle conjugation properties, CP 
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violation, absolute scale of neutrino mass, the mass hierarchy, implications 
for dark matter and large-scale structure, connections with baryogenesis - 
remaining to be resolved. This is one of those wonderful times in physics 
where not only are new discoveries in hand, but we know we have the capacity 
experimentally to seek solutions to even deeper questions. 

2.1 The Standard Solar Model 10 

Observations of stars reveal a wide variety of stellar conditions, with lumi- 
nosities relative to solar spanning a range L ~ 10 -4 to IQPLq and surface 
temperatures T s <~ 2000 — 50000K. The simplest relation one could propose 
between luminosity and T s is that for a blackbody 

L = 47ri? 2 crT 4 => 

(#) 2 #) 4 , (3) 

which suggests that stars of a similar structure might lie along a one-parameter 
path (corresponding to R/ Rq above) in the luminosity (or magnitude) vs. tem- 
perature (or color) plane. In fact, there is a dominant path in the Hertzsprung- 
Russell color-magnitude diagram along which roughly 80% of the stars lie. This 
is the main sequence, those stars supporting themselves by hydrogen burning 
through the pp chain or CNO cycles. 

As one such star, the sun is an important test of our theory of main 
sequence stellar evolution: its properties - age, mass, surface composition, lu- 
minosity, and helioseismology - are by far the most accurately known among 
the stars. The SSM traces the evolution of the Sun over the past 4.6 billion 
years of main sequence burning, thereby predicting the present-day tempera- 
ture and composition profiles of the solar core that govern neutrino production. 
Standard solar models share four basic assumptions: 

• The sun evolves in hydrostatic equilibrium, maintaining a local balance 
between the gravitational force and the pressure gradient. To describe this 
condition in detail, one must specify the equation of state as a function of 
temperature, density, and composition. 

• Energy is transported by radiation and convection. While the solar envelope 
is convective, radiative transport dominates in the core region where ther- 
monuclear reactions take place. The opacity depends sensitively on the solar 
composition, particularly the abundances of heavier elements. 

• Thermonuclear reaction chains generate solar energy. The standard model 
predicts that over 98% of this energy is produced from the pp chain conversion 
of four protons into 4 He (see fig. 1), with proton burning through the CNO 
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cycle contributing the remaining 2%. The Sun is a large but slow reactor: the 
core temperature, T c ~ 1.5 - 10 7 K, results in typical center-of-mass energies for 
reacting particles of ~ 10 keV, much less than the Coulomb barriers inhibiting 
charged particle nuclear reactions. Thus reaction cross sections are small: in 
most cases, as laboratory measurements are only possible at higher energies, 
cross section data must be extrapolated to the solar energies of interest. 
• The model is constrained to produce today's solar radius, mass, and lumi- 
nosity. An important assumption of the standard model is that the Sun was 
highly convective, and therefore uniform in composition, when it first entered 
the main sequence. It is furthermore assumed that the surface abundances of 
metals (nuclei with A > 5) were undisturbed by the subsequent evolution, and 
thus provide a record of the initial solar metallicity. The remaining param- 
eter is the initial 4 Hc/H ratio, which is adjusted until the model reproduces 
the present solar luminosity after 4.6 billion years of evolution. The resulting 
4 Hc/H mass fraction ratio is typically 0.27 ± 0.01, which can be compared 
to the big-bang value of 0.23 ± 0.01. Note that the Sun was formed from 
previously processed material. 

p + p -+ 2 R + e + + v p + p + e -^ 2 H + v 



99.75% 



0.25% 



2 H + p ^ 3 He + 7 



86% 
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3 He + 3 He -^ 4 He + 2p 



3 He + 4 He -^ 7 Be + 7 




7, 



Be + e — > Li + v 
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Figure 1: The solar pp chain. 
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The model that emerges is an evolving Sun. As the core's chemical com- 
position changes, the opacity and core temperature rise, producing a 44% 
luminosity increase since the onset of the main sequence. The temperature 
rise governs the competition between the three cycles of the pp chain: the 
ppl cycle dominates below about 1.6 TO 7 K; the ppll cycle between (1.7-2.3) 
•10 7 K; and the ppIII above 2.4 -10 7 K. The central core temperature of today's 
SSM is about 1.55 -10 7 K. 

The competition between the cycles determines the pattern of neutrino 
fluxes. Thus one consequence of the thermal evolution of our sun is that the 
8 B neutrino flux, the most temperature-dependent component, proves to be 
of relatively recent origin: the predicted flux increases exponentially with a 
doubling period of about 0.9 billion years. 

A final aspect of SSM evolution is the formation of composition gradi- 
ents on nuclear burning timescales. Clearly there is a gradual enrichment of 
the solar core in 4 He, the ashes of the pp chain. Another clement, 3 He, is 
produced and then consumed in the pp chain, eventually reaching some equi- 
librium abundance. The timescale for equilibrium to be established as well as 
the eventual equilibrium abundance are both sharply decreasing functions of 
temperature, and thus increasing functions of the distance from the center of 
the core. Thus a steep 3 He density gradient is established over time. 

The SSM has had some notable successes. From helioseismology the sound 
speed profile c(r) has been very accurately determined for the outer 90% of the 
Sun, and is in excellent agreement with the SSM. Such studies verify important 
predictions of the SSM, such as the depth of the convective zone. However the 
SSM is not a complete model in that it does not explain all features of solar 
structure, such as the depletion of surface Li by two orders of magnitude. This 
is usually attributed to convective processes that operated at some epoch in 
our sun's history, dredging Li to a depth where burning takes place. 

The principal neutrino-producing reactions of the pp chain and CNO cycle 
are summarized in Table 1. The first six reactions produce (3 decay neutrino 
spectra having allowed shapes with endpoints given by E™ ax . Deviations from 
an allowed spectrum occur for 8 B neutrinos because the 8 Be final state is 
a broad resonance. The last two reactions produce line sources of electron 
capture neutrinos, with widths ~ 2 keV characteristic of the temperature of 
the solar core. Measurements of the pp, 7 Be, and 8 B neutrino fluxes will 
determine the relative contributions of the ppl, ppll, and ppIII cycles to solar 
energy generation. As discussed above, and as later illustrations will show more 
clearly, the competition between these cycles is governed in large classes of solar 
models by a single parameter, the central temperature T c . The flux predictions 
of the 1998 calculations of Bahcall, Basu, and Pinsonneault 10 (BP98) and of 
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Table 1: Solar neutrino sources and the flux predictions of the Bahcall/Pinsonneault (BP98) 
and Brun/Turck-Chieze/Morel (BTCM98) SSMs in cm^s" 1 . 



Source 




BP98 


BTCM98 


P + P 


^ 2 H + c+ + v 


0.42 


5.94E10 


5.98E10 


13 N- 


13 C + c+ + v 


1.20 


6.05E8 


4.66E8 


15Q _ 


15 N + c+ + v 


1.73 


5.32E8 


3.97E8 




17 + c+ + v 


1.74 


6.33E6 




8 B -» 


s Bo + c+ + v 


- 15 


5.15E6 


4.82E6 


3 He + p 


-^ 4 Hc + c+ + v 


18.77 


2.10E3 




7 Bc + 


c~ -^ 7 Li + v 


0.86 (90%) 


4.80E9 


4.70E9 






0.38 (10%) 






p + e~ 


+ p -^ 2 H + v 


1.44 


1.39E8 


1.41E8 



Bran, Turck-Chieze and Morel 11 are included in Table 1. 
2.2 Solar Neutrino Detection 12 

Let us start with a brief reminder about low energy neutrino-nucleus interac- 
tions in detectors. Consider the charged current reaction 

u e + {A,Z)^e~ + {A,Z + l) (4) 

Because the momentum transfer to the nucleus is very small for solar neutri- 
nos, it can be neglected in the weak propagator, leading to an effective contact 
current-current interaction. If we begin with the simplest example of a semilep- 
tonic weak process, the decay of a free neutron n rightarrow p+e~ + D e , the 
corresponding transition amplitude is then 

G 

S fi = -JL cos6» c u(p) 7m (1 - g A 75)u(n)u(e) 7 Ai (l - 75)11(1/) (5) 

where Gf is the weak coupling constant measured in muon decay and cos0 c 
gives the amplitude for the weak interaction to connect the u quark to its 
first-generation partner, the d quark. The origin of this effective amplitude 
is the underlying standard model predictions for the elementary quark and 
lepton currents, which are exactly left handed. Experiment shows that the 
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effective coupling of the W boson to the nucleon is governed by 7^(1 — 3.475), 
as noted above, where qa ~ 1.26. The axial coupling is thus shifted from its 
underlying value by the strong interactions responsible for the binding of the 
quarks within the nucleon. 

If an isolated nucleon were the target, one could proceed to calculate the 
cross section from the effective nucleon current given above. The extension 
to nuclear systems traditionally begins with the observation that nucleons in 
the nucleus are rather non-relativistic, v/c — 0.1. The amplitude u(p)7 M (l — 
3A75)w(n) can be expanded in powers of p/M. The leading vector and axial 
operators are readily found to be 



7o 

7 

7075 
775 



1 

p/M -v/c 
a ■ p/M — v/c 
S 



(6) 



Thus it is the time-like part of the vector current and the spacedike part of 
the axial-vector current that survive in the non-relativistic limit." 

If such a non-relativistic reduction is done for our singlc-nucleon current 
one obtains 

Sfi - cos0 c ^t(p)^(n)u(e)7 O (l - <ys)u(i/) 

-0t(p)gA^(n) • u(e)7(l - ts)u(iz)) (7) 

where the <f>'s are now two-component Pauli spinors for the nucleons. The 
above result can be extended to include v e reactions by introducing the isospin 
operators t± where t + | n) = | p) and r_ | p) = | n) , with all other matrix 

"In a nucleus these currents must be corrected for the presence of meson exchange contribu- 
tions. The corrections to the vector charge and axial three-current, which we just pointed 
out survive in the non-relativistic limit, are of order (v/c) 2 ~ 1%. Thus the naive one- 
body currents are a very good approximation to the nuclear currents. In contrast, exchange 
current corrections to the axial charge and vector three-current operators are of order v/c, 
and thus of relative order 1. This difficulty for the vector three-current can be largely cir- 
cumvented, because current conservation as embodied in the generalized Siegert's theorem 
allows one to rewrite important parts of this operator in terms of the vector charge oper- 
ator. In the long-wavelength limit appropriate to solar neutrinos, all terms unconstrained 
by current conservation vanish. In effect, one has replaced a current operator with large 
two-body corrections by a charge operator with only small corrections. In contrast, the 
axial charge operator is significantly altered by exchange currents even for long-wavelength 
processes like /3 decay. Typical axial-charge j3 decay rates are enhanced by ~ 2 because of 
exchange currents. 
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elements being zero. Thus we can generalize our n^p amplitude to n^p by 

0t(p)#n) - ^( N ) T± 0( N ) 

+ (p)CT0(n) -» 0t(N)ff T± 0(N). 

This result easily generalizes to nuclear decay. Given our comments about 
exchange currents, the first step is the replacement 

A 
i=l 

<tt± -> ^<r(i)r±(i). 

i=l 

Plugging S'/j into the standard cross section formula (which involves an average 
over initial and sum over final nuclear spins of the square of the transition 
amplitude) then yields the allowed squared nuclear matrix element 

(l(.fllE T ±«IN)| 2 +5iK/HE^±«IN)| 2 ) • (8) 

1 \ i—\ i=l / 

The Fermi operator is proportional to the isospin raising/lowering opera- 
tor: in the limit of good isospin, which typically is broken at the < 5% level 
for transitions between well-bound nuclear states, the Fermi operator only 
connects states in the same isospin multiplet, that is, states with a common 
spin-spatial structure. If the initial state has isospin (Ti,Mxi), this final state 
has (Ti,Mxi ± 1) for {v, e~) and (D, e + ) reactions, respectively, and is called 
the isospin analog state (IAS). In the limit of good isospin the sum rule for 
this operator in then particularly simple 



2J, 



A A 



The excitation energy of the IAS relative to the parent ground state can be 
estimated accurately from the Coulomb energy difference 

/ i 7987 \ 

Eias ~ 7773 1-293 MeV. (10) 

b VI. 12,4 V3-I-0.78 / 

The angular distribution of the outgoing electron for a pure Fermi (TV, Z) + v — > 
(N — 1, Z + 1) + e~ transition is 1 + (3coa6 e , and thus forward peaked. Here 
/? is the electron velocity. 
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The Gamow- Teller (GT) response is more complicated, as the operator 
connects the ground state to many states in the final nucleus. In general 
we do not have a precise probe of the nuclear GT response apart from weak 
interactions themselves. However a good approximate probe is provided by 
forward-angle (p,n) scattering off nuclei. The (p,n) studies demonstrate that 
the GT strength tends to concentrate in a broad resonance centered at a po- 
sition S = Eqt — Ej as relative to the IAS given by 

^7.0-28.9 ^" MeV. (11) 

Thus while the peak of the GT resonance is substantially above the IAS for 
N ~ Z nuclei, it drops with increasing neutron excess, with (5^0 for Pb. A 
typical value for the full width at half maximum r is ~ 5 MeV. The angular 
distribution of GT (N, Z) + u e -> (N - 1, Z + 1) + e~ reactions is 3 - /3cos6» e , 
corresponding to a gentle peaking in the backward direction. 

The above discussion of allowed responses can be repeated for neutral cur- 
rent processes such as (y,i/). The analog of the Fermi operator contributes 
only to elastic processes, where the standard model nuclear weak charge is 
approximately the neutron number. As this operator does not generate tran- 
sitions, it is not yet of much interest for solar or supernova neutrino detection, 
though there are efforts to develop low-threshold detectors (e.g., cryogenic 
technologies) for recording the modest nuclear recoil energies. The analog of 
the GT response involves 

iM^y)i 2 = ^^i</iii>«^iN}i 2 . (12) 

1 i=l 

The operator appearing in this expression is familiar from magnetic moments 
and magnetic transitions, where the large isovector magnetic moment (fx v ~ 
4.706) often leads to it dominating the orbital and isoscalar spin operators. 

Finally, there is one purely leptonic reaction of great interest, since it is the 
reaction exploited by Kamiokandc and SupcrKamiokande. Electron neutrinos 
can scatter off electrons via both charged and neutral current reactions. The 
cross section calculation is straightforward and will not be repeated here. Two 
features of the result are of importance for our later discussions, however. 
Because of the neutral current contribution, heavy-flavor [y^ and v T ) also 
scatter off electrons, but with a cross section reduced by about a factor of 
seven at low energies. Second, for neutrino energies well above the electron 
rest mass, the scattering is sharply forward peaked. Thus this reaction allows 
one to exploit the position of the Sun in separating the solar neutrino signal 
from a large but isotropic background. 
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As we mentioned earlier, the first experiment performed exploited the re- 
action 

37 Cl(^e~) 37 Ar. 

As the threshold for this reaction is 0.814 MeV, the important neutrino sources 
are the 7 Be and 8 B fluxes. The 7 Be neutrinos excite just the GT transition 
to the ground state, the strength of which is known from the electron capture 
lifetime of 37 Ar. The 8 B neutrinos can excite all bound states in 37 Ar, including 
the dominant transition to the IAS residing at an excitation of 4.99 MeV. The 
strength of excite-state GT transitions can be determined from the (5 decay 
37 Ca(/3 + ) 37 K, which is the isospin mirror reaction to 37 C1(V, e~) 37 Ar. The 
net result is that, for SSM fluxes, 78% of the capture rate should be due 
to 8 B neutrinos, and 15% to 7 Bc neutrinos. The measured capture rate 13 
2.56 ±0.16 ± 0.16 SNU (1 SNU = 10~ 36 capturc/atom/scc) is about 1/3 the 
standard model value. 

Similar radiochemical experiments were done by the SAGE, GALLEX, 
and GNO collaborations using a different target, 71 Ga. The special properties 
of this target include its low threshold and an unusually strong transition to 
the ground state of 71 Ge, leading to a large pp neutrino cross section (see fig. 
2). The experimental capture rates are 66±tt± 3 3 7 2 and 74.1 j^;^;" SNU for 
the SAGE (data through December, 2001) and GALLEX/GNO' (data through 
GNOI) detectors, respectively. The SSM prediction is about 130 SNU 14 . Most 
important, since the pp flux is directly constrained by the solar luminosity in 
all steady-state models, there is a minimum theoretical value for the capture 
rate of 79 SNU, given standard model weak interaction physics. Uncertainties 
in the 71 Ga cross section due to 7 Be neutrino capture to two excited states of 
unknown strength were greatly reduced by direct calibrations of both detectors 
using 51 Cr neutrino sources. 

The water Cerenkov experiments Kamiokande II/III and SuperKamiokande 
viewed solar neutrinos on an event-by-event basis. Solar neutrinos scatter off 
electrons, with the recoiling electrons producing the Cerenkov radiation that 
is then recorded in surrounding photo-tubes. Thresholds are determined by 
background rates; SuperKamiokande operated with triggers as low as 5 MeV. 
The initial experiment, Kamiokande II/III, found a flux of 8 B neutrinos of 
(2.80 ±0.19 ±0.33) • 10 6 /cm 2 s after about a decade of measurement. Its much 
larger successor SuperKamiokande, with a 22.5 kiloton fiducial volume, yielded 
the result (2.35 ± 0.02 ± 0.08) • 10 6 /cm 2 s after 1496 days of measurements 15 , 
corresponding to 0.465 of the SSM 8 B neutrino flux. 

Results from the Sudbury Neutrino Observatory (SNO) will be discussed 
later in the lectures. 



11 



'Ge + n™ 6 



500 


3/2" 


175 


5/2" 





1/2" 



71 



Ge 



o 



71 



Ga 



pp 7 Be 51 Cr 8 B 
neutrino sources 



Figure 2: Level scheme for Gc showing the excited states that contribute to absorption of 
pp, 7 Be, 51 Cr, and 8 B neutrinos. 



2.3 The Argument for New Particle Physics 

The pattern of solar neutrino fluxes that has emerged from these experiments 
is 

0(pp)^O.9 SSM (pp) 
</>( 7 Be) ~ 

(t)( 8 B) ~ O.43 SSM ( 8 B). (13) 

A reduced 8 B neutrino flux can be produced by lowering the central temper- 
ature of the sun somewhat, as </>( 8 B)^ T c 18 . However, such an adjustment, 
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either by varying the parameters of the SSM or by adopting some nonstan- 
dard physics, tends to push the 0( 7 Be)/(/>( 8 B) ratio to higher values rather 
than the low one of eq. (13), 

„ (14) 

Thus the observations seem difficult to reconcile with plausible solar model 
variations: one observable, <j>{ 8 B), requires a cooler core while a second, the 
ratio </>( 7 Be)/0( 8 B), requires a hotter one. 

How robust is this apparent contradiction? In the last decade, as evidence 
mounted that the solar neutrino problem was a profound one, a key issue was 
SSM uncertainties. Inputs into the SSM - pp chain nuclear cross sections, 
solar parameters like the age, luminosity, and composition, and atomic physics 
quantities such as opacities and screening corrections - all have measurement 
and theory errors. It was the care with which these uncertainties were assessed 
that convinced the community that the solar neutrino problem was a serious 
one. 

No issue received more scrutiny that the nuclear cross sections. The pp 
chain involves a series of non-resonant charged-particle reactions occurring at 
ccntcr-of-mass energies that are well below the height of the inhibiting Coulomb 
barriers. As the resulting small cross sections generally preclude laboratory 
measurements at the relevant energies, one must extrapolate higher energy 
measurements to threshold to obtain solar cross sections. This extrapolation 
is usually discussed in terms of the astrophysical S-factor 

CT (£) = ^exp(-27n?) (15) 

where r\ = Zl ^ 2a ; with a the fine structure constant and f3 = v/c the relative 
velocity of the colliding particles. This parameterization removes the gross 
Coulomb effects associated with the s-wave interactions of charged, point-like 
particles. The remaining energy dependence of S(E) is gentle and can be ex- 
pressed as a low-order polynomial in E. Usually the variation of S(E) with E is 
taken from a direct reaction model and then used to extrapolate higher energy 
measurements to threshold. The model accounts for finite nuclear size effects, 
strong interaction effects, contributions from other partial waves, etc. As labo- 
ratory measurements are made with atomic nuclei while conditions in the solar 
core guarantee the complete ionization of light nuclei, additional corrections 
must be made to account for the different electronic screening environments. 

Thus a great deal of effort was invested in laboratory measurements, in 
the theory required to extrapolate those measurements to the Gamow peak 



13 



energies where solar reactions take place, and in assessing the resulting cross 
section uncertainties. In addition, other SSM uncertainties were evaluated. 
One qualitative illustration of the results is given by fig. 3, which summa- 
rizes a Monte Carlo study of SSM input parameter uncertainties. Five key 
input parameters, the primordial heavy-element-to- hydrogen ratio Z/X and 
S(0) for the p-p, 3 He- 3 He, 3 He- 4 He, and p- 7 Be were varied according to their 
assigned uncertainties, assuming for each parameter a normal distribution with 
the appropriate mean and standard deviation. (These five were the parame- 
ters assigned the largest uncertainties.) Smaller uncertainties from radiative 
opacities, the solar luminosity, and the solar age were folded into the results 
of the model calculations perturbatively 16 > 17 . 

The resulting pattern of Be and 8 B flux predictions produces an elongated 
error ellipse, showing that v flux changes are strongly correlated even when 
a large set of distinct uncertainties are explored. Those variations producing 
</>( 8 B) below O.80 SSM ( 8 B) tend to produce a reduced 0( 7 Be), but the reduction 
is always less than 0.8. Thus a greatly reduced (j)( 7 Be) cannot be achieved 
within the uncertainties assigned to parameters in the SSM. 

A similar exploration, but including parameter variations very far from 
their preferred values, was carried out by Castellani et al. 18 , who displayed 
their results as a function of the resulting core temperature T c . The pattern 
that emerges is striking (see fig. 4): parameter variations producing the same 
value of T c produce remarkably similar fluxes. Thus T c provides an excellent 
one-parameter description of standard model perturbations. Figure 4 also 
illustrates the difficulty of producing a low ratio of 0( 7 Be)/</>( 8 B) when T c is 
reduced. 

The Monte Carlo parameter variations of fig. 3 were constrained to repro- 
duce the solar luminosity. Those variations show a similar strong correlation 
with T c 

^(pp)aT- 1 - 2 0( 7 Bc)cxT 8 <f>( 8 B) ex T c 18 . (16) 

Figures 3 and 4 offer a strong argument that reasonable variations in the pa- 
rameters of the SSM, or nonstandard changes in quantities like the metallicity, 
opacities, or solar age, cannot produce the pattern of fluxes deduced from 
experiment (eq. (13)). This would seem to limit possible solutions to errors 
either in the underlying physics of the SSM or in our understanding of neutrino 
properties. 

The Castellani et al. explorations belong to a larger class of proposed non- 
standard solar models where either very large parameter changes (in nuclear 
cross sections, opacities, etc.) or new physics (e.g., mixing in the solar core) are 
hypothesized. Once results from CI, SAGE/GALLEX, and Kamiokande were 
available, it became possible to argue that no such nonstandard solar model 
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Figure 3: SSM 7 Be and S B flux predictions. The dots represent the results of SSM calcu- 
lations where the input parameters were varied according to their assigned uncertainties, as 
described in the text. The 90% and 99% confidence level error ellipses are shown. 



can solve the solar neutrino problem: if one assumes undistorted neutrino spec- 
tra, no combination of pp, 7 Be, and 8 B neutrino fluxes fits the experimental 
results well 19 . In fact, in an unconstrained fit, the required 7 Be flux is un- 
physical, negative by almost 3o\ Thus, barring some unfortunate experimental 
error, new particle physics seemed to be the indicated solution. This conclu- 
sion was reinforced by the successes the SSM had in reproducing the new data 
on helioseismology 20 . Suggested particle physics solutions include neutrino os- 
cillations, neutrino decay, neutrino magnetic moments, and weakly interacting 
massive particles. Among these, the Mikheyev-Smirnov-Wolfenstein effect — 
neutrino oscillations enhanced by matter interactions — was widely regarded 
as the most plausible. 
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Figure 4: The responses of the pp, 7 Be, and 8 B neutrino fluxes to the indicated variations in 
solar model input parameters, displayed as a function of the resulting central temperature 
T c . From Castcllani et al. 



2.4 Neutrino Oscillations 

One odd feature of particle physics is that neutrinos, which are not required 
by any symmetry to be massless, nevertheless must be much lighter than any 
of the other known fermions. For instance, the current limit on the v e mass 
is < 2.2 eV. The standard model requires neutrinos to be massless, but the 
reasons are not fundamental. Dirac mass terms mn, analogous to the mass 
terms for other fermions, cannot be constructed because the model contains 
no right-handed neutrino fields. Neutrinos can also have Majorana mass terms 

v c L m,LVL and ^hirz/r (17) 
16 



where the subscripts L and R denote left- and right-handed projections of 
the neutrino field v 1 and the superscript c denotes charge conjugation. The 
first term above is constructed from left-handed fields, but can only arise as 
a nonrenormalizable effective interaction when one is constrained to generate 
mi with the doublet scalar field of the standard model. The second term is 
absent from the standard model because there are no right-handed neutrino 
fields. 

None of these standard model arguments carries over to the more general, 
unified theories that theorists believe will supplant the standard model. In the 
enlarged multiplets of extended models it is natural to characterize the fermions 
of a single family, e.g., v e , c, u, d, by the same (Dirac) mass scale mo- Indeed 
the charged members of this multiplet all have comparable masses ~ 1 MeV. 
The neutrino, however, is much lighter; but the neutrino, which carries no 
charge, is the only standard model fermion that can have both Majorana and 
Dirac mass terms. Thus it is natural to explain small neutrino masses as a 
consequence of Majorana masses. In the seesaw mechanism? 1 , 

M v ~(\ mD ). (18) 

2 

Diagonalization of the mass matrix produces one light neutrino, mught ~ 
and one unobservably heavy, Whoavy ~ TOr. The factor (m_o/m^) is the needed 
small parameter that accounts for the distinct scale of neutrino masses. The 
masses for the f e ,^, and v T are then related to the squares of the corre- 
sponding quark masses m u , m c , and m t . Taking tur <~ 10 16 GcV for the 
right-handed Majorana mass, a typical grand unification scale for models built 
on groups like SO(10), the seesaw mechanism gives the crude relation 

m Ve : m Vlx : m„ T <-> 2 • 1CT 12 : 2 • 1(T 7 : 3 • l(T 3 eV. (19) 

The fact that solar neutrino experiments can probe small neutrino masses, and 
thus provide insight into possible new mass scales vhr that are far beyond the 
reach of direct accelerator measurements, has been an important theme of the 
field. 

One of the most interesting possibilities for solving the solar neutrino prob- 
lem has to do with neutrino masses. For simplicity we will discuss just two 
neutrinos. If a neutrino has a mass to, we mean that as it propagates through 
free space, its energy and momentum are related in the usual way for this mass. 
Thus if we have two neutrinos, we can label those neutrinos according to the 
eigenstates of the free Hamiltonian, that is, as mass eigenstates. 

But neutrinos are produced by the weak interaction. In this case, we have 
another set of eigenstates, the flavor eigenstates. We can define a v e as the 
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neutrino that accompanies the positron in f3 decay. Likewise we label by 
the neutrino produced in muon decay. 

Now the question: are the eigenstates of the free Hamiltonian and of the 
weak interaction Hamiltonian identical? Most likely the answer is no: we know 
this is the case with the quarks, since the different families (the analog of the 
mass eigenstates) do interact through the weak interaction. That is, the up 
quark decays not only to the down quark, but also occasionally to the strange 
quark. (This is why we had a cos 9 C in our weak interaction amplitude: the 
amplitude for u —* s is proportional to sin# c .) Thus we suspect that the weak 
interaction and mass eigenstates, while spanning the same two-neutrino space, 
are not coincident: the mass eigenstates \i/i) and \v 2 ) (with masses mi and 
m-i) are related to the weak interaction eigenstates by 

\y e ) = cos6 v \v\) +sm.6 v \wi) 

Wfj.) —- sin 9 v \vi) + cos 6 v \v 2 ) (20) 

where 9 V is the (vacuum) mixing angle. 

An immediate consequence is that a state produced as a \v e ) or a 
at some time t — for example, a neutrino produced in [5 decay — does not 
remain a pure flavor eigenstate as it propagates away from the source. The 
different mass eigenstates comprising the neutrino will accumulate different 
phases as they propagate downstream, a phenomenon known as vacuum os- 
cillations (vacuum because the experiment is done in free space). To see the 
effect, suppose the neutrino produced in a (3 decay is a momentum eigenstate. 
At time t=0 

\u(t = 0)) = \v e ) =cos0„|i/i) +sm6 v \v 2 ). (21) 
Each eigenstate subsequently propagates with a phase 

gi(k-x-ujt) ^i(k-x — ^/m^+k 2 t) (22) 

But if the neutrino mass is small compared to the neutrino momentum/energy, 
one can write 

V / mf + fc^fc(l + g). (23) 

Thus we conclude 

\u{t)) — ^(k-Z-kt-iml+m^t/Ak) 

x[cos9 v \ Vl )e tSm2t ' Ak + sin9 v \v 2 )e- tSm2t ' Ak ]. (24) 

We see there is a common average phase (which has no physical consequence) 
as well as a beat phase that depends on 

5m 2 = m\ — m 1 . (25) 
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Now it is a simple matter to calculate the probability that our neutrino state 
remains a \u e ) at time t 
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P Ve {t) = \{y e Ht))\ 

= 1 - sin 2 29 v sin 2 (^j^j 1 - \ sin 2 29 v (26) 

where the limit on the right is appropriate for large t. (When one properly de- 
scribes the neutrino state as a wave packet, the large-distance behavior follows 
from the eventual separation of the mass eigenstates.) Now E <~ k, where E 
is the neutrino energy, by our assumption that the neutrino masses are small 
compared to k. We can reinsert the implicit constants to write the probability 
in terms of the distance x of the neutrino from its source, 

P,(x) = l-sin 2 2^sin 2 (^^). (27) 

If the oscillation length 

Lo = -hh (28) 

is comparable to or shorter than one astronomical unit, a reduction in the solar 
v e flux would be expected in terrestrial neutrino oscillations. 

The suggestion that the solar neutrino problem could be explained by 
neutrino oscillations was first made by Pontecorvo in 1958, who pointed out 
the analogy with Kq <-» Kq oscillations. From the point of view of particle 
physics, the sun is a marvelous neutrino source. The neutrinos travel a long 
distance and have low energies (~ 1 MeV), implying a sensitivity down to 

Sm 2 > lCr 12 eV 2 . (29) 

In the seesaw mechanism, Sm 2 ~ m 2 ., so neutrino masses as low as m 2 ~ 1CP 6 
eV could be probed. 

From the expressions above one expects vacuum oscillations to affect all 
neutrino species equally, if the oscillation length is small compared to an as- 
tronomical unit. This is somewhat in conflict with the solar neutrino data, as 
we have argued that the 7 Be neutrino flux is quite suppressed. Furthermore, 
there is a weak theoretical prejudice that 6 V should be small, like the Cabibbo 
angle. The first objection, however, can be circumvented in the case of "just 
so" oscillations where the oscillation length is comparable to one astronomical 
unit. In this case the oscillation probability becomes sharply energy dependent, 
and one can choose Sm 2 to preferentially suppress one component (e.g., the 
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monochromatic 7 Be neutrinos), though the requirement for large mixing an- 
gles remains. This "just so" vacuum scenario is one that received considerable 
attention. 

Below we will find that oscillations in matter can lead to nearly total flavor 
conversion even if the mixing angle is small. In preparation for this we first 
present the results above in a slightly more general way. The analog of eq. 
(24) for an initial muon neutrino (\v(t = 0)) = 1^)) is 

\v(t)) = ^(k-Z-kt-iml+mDt I '4k) 

x[-sin0>i)e l5m2t/4 ' £ +cosd v \is 2 }e- lSm2t/4k } (30) 

Now if we compare Eqs. (24) and (30) we see that they are special cases of a 
more general problem. Suppose we write our initial neutrino wave function as 

\v{t = 0)) - a e (t = 0)K) + a^t = 0)|^). (31) 

Then Eqs. (24) and (30) tell us that the subsequent propagation is described 
by changes in a e (x) and a^x) according to (this takes a bit of algebra) 

'— ( ae \ — — ( ~3m 2 cos2# w <5m 2 sin 20 v \ ( a e \ ,„„, 

l dx V a M/ 4^ V Sm 2 sm26 v 6m 2 cos 29 v ) \ a M ) ' ' ' 

Note that the common phase has been ignored: it can be absorbed into the 
overall phase of the coefficients a e and a M , and thus has no consequence. Also, 
we have equated x = t, that is, set c = 1. 



2.5 The Mikheyev-Smirnov-Wolfenstein Mechanism 

The view of neutrino oscillations changed when Mikheyev and Smirnov 22 
showed in 1985 that the density dependence of the neutrino effective mass, 
a phenomenon first discussed by Wolfcnstcin in 1978, could greatly enhance 
oscillation probabilities: a v e is adiabatically transformed into a as it tra- 
verses a critical density within the sun. It became clear that the sun was not 
only an excellent neutrino source, but also a natural regenerator for cleverly 
enhancing the effects of flavor mixing. 

While the original work of Mikheyev and Smirnov was numerical, their 
phenomenon was soon understood analytically as a level-crossing problem. If 
one writes the neutrino wave function in matter as in eq. (31), the evolution 
of a e (x) and a t _ l (x) is governed by 

.d_fa e \ _ J_ / 2EV2G F p(x) - Sm 2 cos 20 v 5m 2 sin 28 v \ / a e \ 
1 dx V «m / ~ AE V s m 2 sin 20 v - 2E^/2G F p{x) + dm 2 cos 26 v ) \ ) 

(33) 



20 



where Gf is the weak coupling constant and p(x) the solar electron density. 
If p(x) = 0, this is exactly our previous result and can be trivially integrated 
to give the vacuum oscillation solutions given above. The new contribution 
to the diagonal elements, 2E\/2Gfp{x), represents the effective contribution 
to the M 2 matrix that arises from neutrino-electron scattering. The indices 
of refraction of electron and muon neutrinos differ because the former scatter 
by charged and neutral currents, while the latter have only neutral current 
interactions. The difference in the forward scattering amplitudes determines 
the density-dependent splitting of the diagonal elements of the matter equation, 
the generalization of cq. (32). 

It is helpful to rewrite this equation in a basis consisting of the light and 
heavy local mass eigenstates (i.e., the states that diagonalize the right-hand 
side of eq. (33)), 

\vl{x)) = cos9(x)\v e ) - sm0(x)\v^) 

\v H {x)) = sva.e{x)\v e ) +cos6{x)\v IJ ). (34) 

The local mixing angle is defined by 

sin 29 v 



sin2#(x) 



sjx 2 (x)+sm 2 29 v 



cos2fl(x) = X ^ = (35) 



sjx 2 (x) +sin 2 26v 



where X(x) = 2\[2Gfp(x)E / '5m 2 — cos2# w . Thus 9(x) ranges from 9 V to tt/2 
as the density p(x) goes from to oo. 
If we define 

\y{x)) = a H (x)\v H (x)) + a L (x)\v L (x)), (36) 
the neutrino propagation can be rewritten in terms of the local mass eigenstates 

i *(*") = ( A (f ia ^)( aH ) (37) 
dx \a L I \-ia(x) -\{x) I \a L I 



with the splitting of the local mass eigenstates determined by 

dm 2 
~2E 

and with mixing of these eigenstates governed by the density gradient 



2A(x) = ^ V /X2(x)+Sin 2 20 v (38) 



_ f M V2GF±p(x)sin29 v 
a{X) -{sm 2 ) X 2 (x)+sin 2 29 v ' ^ 
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The results above are quite interesting: the local mass eigenstates diagonalize 
the matrix if the density is constant. In such a limit, the problem is no more 
complicated than our original vacuum oscillation case, although our mixing an- 
gle is changed because of the matter effects. But if the density is not constant, 
the mass eigenstates evolve as the density changes. This is the crux of the 
MSW effect. Note that the splitting achieves its minimum value, sin 26 v , 
at a critical density p c = p{x c ) 

2\p2EG F p c = 5m 2 cos 26 v (40) 

that defines the point where the diagonal elements of the original flavor matrix 
cross. 

Our local-mass-eigenstate form of the propagation equation can be trivially 
integrated if the splitting of the diagonal elements is large compared to the off- 
diagonal elements (see eq. (37)), 



X(x) 



a(x) 



sin 2 29 v 6m 2 1 [X{xf + sin 2 29 v f/ 2 ^ 
cos26» v ~2E |J_£W£)| sin 3 26 v > ' 



a condition that becomes particularly stringent near the crossing point, 



7c = j(x c ) 



sin 2 26 v 8m 2 1 
cos 26 v ~2E 



J_ dp(x) I 
p c dx \ x — x c 



> 1. (42) 



That is, adiabaticity depends on the density scale height at the crossing point. 
The resulting adiabatic electron neutrino survival probability 23 , valid when 
7c > 1, is 

P* diab = \ + \ cos 26 v cos 26i (43) 

where 0, = 6{xi) is the local mixing angle at the density where the neutrino 
was produced. 

The physical picture behind this derivation is illustrated in Figure 5. One 
makes the usual assumption that, in vacuum, the v e is almost identical to the 
light mass eigenstate, z^l(O), i.e., m\ < m-i and cosf^ ~ 1. But as the density 
increases, the matter effects make the v e heavier than the f M , with v e — ► vh{x) 
as p(x) becomes large. The special property of the Sun is that it produces u e s 
at high density that then propagate to the vacuum where they are measured. 
The adiabatic approximation tells us that if initially v e <~ vh{x), the neutrino 
will remain on the heavy mass trajectory provided the density changes slowly. 
That is, if the solar density gradient is sufficiently gentle, the neutrino will 
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Figure 5: Schematic illustration of the MSW crossing. The dashed lines correspond to the 
electron-electron and muon-muon diagonal elements of the matrix in the flavor basis. 
Their intersection defines the level-crossing density p c . The solid lines are the trajectories 
of the light and heavy local mass eigenstates. If the electron neutrino is produced at high 
density and propagates adiabatically, it will follow the heavy-mass trajectory, emerging from 
the sun as a v^. 



emerge from the sun as the heavy vacuum eigenstate, ~ v^. This guarantees 
nearly complete conversion of u e s into v^s, producing a flux that cannot be 
detected by the Homestake or SAGE/GALLEX detectors. 

But this does not explain the curious pattern of partial flux suppressions 
coming from the various solar neutrino experiments. The key to this is the 
behavior when j c < 1. Our expression for 7(2;) shows that the critical region for 
non-adiabatic behavior occurs in a narrow region (for small 6 V ) surrounding the 
crossing point, and that this behavior is controlled by the density scale height. 
This suggests an analytic strategy for handling non-adiabatic crossings: one 
can replace the true solar density by a simpler (integrablc!) two-parameter 
form that is constrained to reproduce the true density and its derivative at 
the crossing point x c . Two convenient choices are the linear (p(x) = a + 
bx) and exponential (p(x) — ae~ bx ) profiles. As the density derivative at x c 
governs the non-adiabatic behavior, this procedure should provide an accurate 
description of the hopping probability between the local mass eigenstates when 
the neutrino traverses the crossing point. The initial and ending points Xi and 
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Xf for the artificial profile are then chosen so that p{x{) is the density where 
the neutrino was produced in the solar core and p(xf) = (the solar surface), 
as illustrated in in Figure 6. Since the adiabatic result (P^f lah ) depends only 
on the local mixing angles at these points, this choice builds in that limit. But 
our original flavor-basis equation can then be integrated exactly for linear and 
exponential profiles, with the results given in terms of parabolic cylinder and 
Whittaker functions, respectively. 
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Figure 6: The top figure illustrates, for one choice of sin 2 29 and Sm 2 , that the region of 
non-adiabatic propagation (solid line) is usually confined to a narrow region around the 
crossing point r c . In the lower figure, the solid lines represent the solar density and a 
linear approximation to that density that has the correct initial and final values, as well 
as the correct density and density derivative at r c . Thus the linear profile is a very good 
approximation to the sun in the vicinity of the crossing point. The MSW equations can 
be solved analytically for this wedge. By extending the wedge to ±oo (dotted lines) and 
assuming adiabatic propagation in these regions of unphysical density, one obtains the simple 
Landau-Zcncr result discussed in the text. 

That result can be simplified further by observing that the non-adiabatic 
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region is generally confined to a narrow region around x c , away from the end- 
points Xi and Xf. We can then extend the artificial profile to x = ±00, as 
illustrated by the dashed lines in Figure 6. As the neutrino propagates adi- 
abatically in the unphysical region x < x i: the exact solution in the physical 
region can be recovered by choosing the initial boundary conditions 

. . , . . -i \ Xl X(x)dx 

a,L{— 00) = — a^(— 00) = cosflje J -°° 

an{— 00) = a e (— 00) = smOie x ^ dx _ (44) 

That is, \v{— 00)) will then adiabatically evolve to \v(xi)) — \u e ) as x goes from 
—00 to Xi. The unphysical region x > Xf can be handled similarly. 

With some algebra a simple generalization of the adiabatic result emerges 
that is valid for all Sm 2 /E and 9 V 

P -e = l + l cos 29 v cos 20j(l - 2P hop ) (45) 

where Phop is the Landau-Zener probability of hopping from the heavy mass 
trajectory to the light trajectory on traversing the crossing point. For the 
linear approximation to the density 24 > 25 , 

= e~* Jc/2 - (46) 

As it must by our construction, P„ c reduces to p^" ab for j c 1. When the 
crossing becomes non-adiabatic (e.g., "f c <§; 1 ), the hopping probability goes 
to 1, allowing the neutrino to exit the sun on the light mass trajectory as a i/ e , 
i.e., no conversion occurs. 

Thus there are two conditions for strong conversion of solar neutrinos: 
there must be a level crossing (that is, the solar core density must be sufficient 
to render v e ~ vn(xi) when it is first produced) and the crossing must be 
adiabatic. The first condition requires that Sm 2 /E not be too large, and the 
second 7 C > 1. The combination of these two constraints, illustrated in fig. 

7, defines a triangle of interesting parameters in the sin 29 v plane, 

as Mikheyev and Smirnov found by numerically integration. A remarkable 
feature of this triangle is that strong v e — > conversion can occur for very 
small mixing angles (sin 2 29 ~ 1CT 3 ), unlike the vacuum case. 

One can envision superimposing on fig. 7 the spectrum of solar neutri- 
nos, plotted as a function of ^g- for some choice of 8m 2 . Since Davis sees 
some solar neutrinos, the solutions must correspond to the boundaries of the 
triangle in fig. 7. The horizontal boundary indicates the maximum for 
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Figure 7: MSW conversion for a neutrino produced at the sun's center. The upper shaded 
region indicates those 8m? /E where the vacuum mass splitting is too great to be overcome 
by the solar density. Thus no level crossing occurs. The lower shaded region defines the 
region where the level crossing is non-adiabatic (7 C less than unity). The unshaded region 
corresponds to adiabatic level crossings where strong u e — > will occur. 



which the sun's central density is sufficient to cause a level crossing. If a spec- 
trum properly straddles this boundary we obtain a result consistent with the 
Homestakc experiment in which low energy neutrinos (large 1/E) lie above 
the level-crossing boundary (and thus remain f e 's), but the high-energy neu- 
trinos (small 1/E) fall within the unshaded region where strong conversion 
takes place. Thus such a solution would mimic nonstandard solar models in 
that only the 8 B neutrino flux would be strongly suppressed. The diagonal 
boundary separates the adiabatic and non-adiabatic regions. If the spectrum 
straddles this boundary we obtain a second solution in which low energy neu- 
trinos lie within the conversion region, but the high-energy neutrinos (small 
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1/E) lie below the conversion region and are characterized by 7 <C 1 at the 
crossing density. (Of course, the boundary is not a sharp one, but is charac- 
terized by the Landau-Zener exponential). Such a non-adiabatic solution is 
quite distinctive since the flux of pp neutrinos, which is strongly constrained 
in the standard solar model and in any steady-state nonstandard model by 
the solar luminosity, would now be sharply reduced. Finally, one can imag- 
ine "hybrid" solutions where the spectrum straddles both the level-crossing 
(horizontal) boundary and the adiabaticity (diagonal) boundary for small 9, 
thereby reducing the 7 Be neutrino flux more than either the pp or 8 B fluxes. 




Figure 8: Illustration of the SMA, LMA, and LOW solutions arising out of fits to the event 
rates for CI, Kamiokande, and SAGE/GALLEX. Figure provided by K. Heeger. 
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What are the results of a careful search for MSW solutions satisfying 
the Homestake, Kamiokande, and SAGE/GALLEX constraints? This was 
explored in detail by several groups (see fig. 8). One solution, corresponding 
to a region surrounding 8m 2 <~ 6 • 10~ 6 eV 2 and sin 2 20 v ~ 6 • 10~ 3 , is the hybrid 
case described above. It is commonly called the small-mixing-angle solution 
(SMA). A second, large-angle solution (LMA) corresponds to Sm 2 <~ 10~ 5 eV 2 
and sin 2 26 v ^0.6. These solutions can be distinguished by their characteristic 
distortions of the solar neutrino spectrum. The survival probabilities P^ ISW (E) 
for the small- and large-angle parameters given above are shown as a function 
of E in fig. 9. 



Figure 9: MSW survival probabilities P(E„) for typical small angle and large angle solutions. 

The MSW mechanism provides a natural explanation for the pattern of 
observed solar neutrino fluxes. While it requires profound new physics, both 
massive neutrinos and neutrino mixing are expected in extended models. 

2.6 SuperKamiokande, SNO, and the Neutrino Mixing Matrix 

Over the past five years the Cl/SAGE/GALLEX/Kamiokandc hints of new 
neutrino physics have been spectacularly confirmed by SuperKamiokande and 
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SNO, and a program to further constrain the neutrino mass matrix with ter- 
restrial neutrino experiments has been inaugurated by KamLAND and K2K. 

SupcrKamiokande and Sudbury Neutrino Observatory (SNO) detectors 
are real-time counting detectors, in contrast to the radiochemical detectors 
such as Homestake, GALLEX/GNO, and SAGE, which can only determine a 
time- and energy-integral of the flux. Both SupcrKamiokande and SNO can 
detect neutrinos through elastic scattering 

Vx + e~ — > v x + e~ . (47) 

The electrons coming from this reaction are confined to a forward cone, with 
most of the observed spread in angles around the forward peak coming from 
the limited angular resolution of the detector. To an excellent approxima- 
tion the recorded events (including the effects of resolution) are contained in 
the forward hemisphere. Thus experimentalist can use the isotropic distribu- 
tion of events in the backward hemisphere (defined around the vector pointing 
from the sun) to determine the background, then substract in the forward 
hemisphere to obtain the events coming from solar neutrinos. The neutrino 
energy is difficult to reconstruct because the initial neutrino momentum is 
shared by the scattered neutrino and electron. Nevertheless, for MSW solu- 
tions like the SMA where there is a significant distortion of the v e spectrum, 
there is a sufficient residual distortion of the electron spectrum to signal new 
physics. The SuperKamiokandc collaboration carefully calibrated the detector 
with electrons from a linac so that small spectral distortions could be reliably 
extracted. 

In addition to the reaction eq. (47) SNO can detect neutrinos by two 
additional reactions, one via the charged current 

v e +d^ p + p + e", (48) 

and the second via neutral current scattering 

^x(^x) + d^^ x (F x )+p + n. (49) 

The neutrons produced in eq. (49) can be detected either by (n,7) on the 
heavy water or on a salt introduced to enhance the capture, or by using 3 Hc 
proportional counters. The electrons coming from the reaction (48) are quite 
hard, with energies not too different from <~ E v — f .44 MeV, and with a angular 
distribution approximately that of a pure GT transition in the relativistic limit, 
(1 — cos6* e /3) with respect to the incident neutrino. This backward peaking 
contrasts nicely with the forward-peaked elastic scattering signal. 
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Figure 10: The dashed line shows the spectrum distortion at SNO for the small-angle MSW 
solution (Sm 2 ~ 5 X 10~ 6 eV 2 and sin 20 ~ 0.01). The solid line is the spectrum without 
MSW oscillations, normalized to the same total rate as with MSW oscillations. 



The hardness of the charged current reaction (eq. (48)) makes it an effec- 
tive tool for finding spectral distortions induced by the MSW mechanism. The 
SNO spectral distortion expect for the SMA solution (Sm 2 ~ 5x 1CP 6 eV 2 and 
sin 26 <~ 0.01) is shown in fig. 10. 

If there are no oscillations into sterile states (new neutrino states lacking 
the usual standard model weak interactions), the neutral current reaction (eq. 
(49)) measures the total SSM flux, independent of flavor. One other reaction 
of potential interest 

i> e + d^n + n + e + (50) 

produces a two-neutron coincidence. Electron antineutrinos can arise from 
spin-flavor oscillations 26 ' 27 in the sun and, of course, from supernovae. 
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As the solar neutrino problem deepened with the new measurements from 
SAGE/GALLEX and Kamiokande, a second neutrino anomaly also began to 
draw attention. Several early underground detectors (1MB, Kamiokande, and 
later Soudan II) found a deficit of muons produced by atmospheric neutrinos. 
Atmospheric neutrinos arise from the decay of secondary pions, kaons, and 
muons produced by the collisions of primary cosmic rays with the oxygen and 
nitrogen nuclei in the upper atmosphere. For energies less than 1 GeV all the 
secondaries decay : 

M ± -» e ± + v e (V e ) +V ll {v tl ) 

Consequently one expects the ratio 

r = (v e +V e )/{v^ +77^) (52) 

to be approximately 0.5 in this energy range. Detailed Monte Carlo calcu- 
lations 28 , including the effects of muon polarization, give r ~ 0.45. As one 
is evaluating a ratio of similarly calculated processes, r does not require an 
absolute flux calculation and is thus relative free of systematic uncertainties. 
Different groups estimating this ratio, even though they start with neutrino 
fluxes which can differ in magnitude by up to 25%, all agree within a few 
percent 29 . As the shower energy increases more muons survive due to time 
dilation. Hence one expects the ratio r to decrease as the energy increases. 
The ratio (observed to predicted) of ratios 

_ iyjj_ ^e)data 
iy^l jv e ) MontoCarlo 

studied by the experimentalists should then be unity, in the absence of oscil- 
lations. 

The first "smoking gun" for neutrino oscillations came from the detailed 
atmospheric v results of SupcrKamiokandc. The initial announcement of neu- 
trino oscillations, made in 1998, is now supported by 1489 days of data. The 
experimenters found 

R = 0.638 ± 0.016(stat) ± 0.050(syst) (54) 
for sub-GeV events which were fully contained in the detector and 

R = 0.658tS(stat) ± 0.078(syst) (55) 
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Figure 11: The ratio of fully contained events measured at SupcrKamiokande versus recon- 
structed L/E v . The dashed lines show the expected shape for <-> v T oscillations with 
6m 2 = 2.2 • 1(T 3 eV 2 and sin 2 29 = 1. 



for fully- and partially-contained multi-GeV events. The large deviation from 
R — 1 indicates neutrino oscillations. Much more dramatic evidence for oscil- 
lations comes from measuring R as a function of the zenith angle, 6, between 
the vertical and incident neutrino. A down-going neutrino (0 ~ 0°) trav- 
els through the atmosphere above the detector (a distance of about 20 km), 
whereas an up-going neutrino (9 <~ 180°) has traveled through the entire Earth 
(a distance of about 13000 km). Hence a measurement of the flux as a function 
of the zenith angle yields information about neutrino survival probabilities as 
a function of the distance traveled. 

The SuperKamiokande collaboration measured the zenith angle depen- 
dence not only of R, but also of the electron and muon neutrino fluxes sepa- 
rately 9 . This information is shown in Fig. 11, where the data is plotted as 
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a function of the reconstructed L/E v instead of the zenith angle. The data 
exhibit a zenith-angle (distance) dependent deficit of muon neutrinos, but not 
of electron neutrinos, a behavior consistent with — > v T oscillations. This 
interpretation is consistent with the deficits of up-going muons measured with 
the Kamiokande 30 and MACRO 31 detectors. These muons are produced by 
the very high energy up-going muon neutrinos in the rock surrounding the 
detectors. 

One remarkable feature of the SuperKamiokande atmospheric neutrino 
results is the deduced mixing angle. Designating the two mass eigenstates 
participating in the mixing as 2 and 3, the best fit sin 2 2# 2 3 is 1, with sin 2 #23 > 
0.92 at 90% c.l. As quark mixing angles are small, a v mixing angle near 
maximal (#23 ~ 45°) was not anticipated. The best-fit 6m 2 3 is <~ 2.5 • 10~ 3 
cV 2 . 

This year an equally spectacularly resolution of the solar neutrino problem 
was obtained by SNO. SNO's proof of oscillations was definitive: the heavy- 
flavor neutrinos produced as a result of solar v e oscillations were seen directly, 
by comparing the charge current and neutral current results. 

SNO's initial results arc shown in fig. 12. The detector's three v channels 
- charge and neutral current scattering off deuterium and ^-electron clastic 
scattering - combine to show that approximately two-thirds of the solar flux 
is in heavy flavors. The fluxes deduced, assuming a standard 8 B v spectrum 
shape, are 8 

<t> Ve = (1.76 ± Om(stat) ± 0.09(syst)) ■ 10 6 /cm 2 s 
i, e „ s = (3.41 ± 0A5(stat)t a o it(syst)) • 10 6 /cm 2 s. (56) 

The 5.3 er difference is the significance of the oscillation proof. The total active 
flux, measured above the neutron breakup threshold for deuterium of 2.2 MeV, 
is 

0„ = 5.09±g;t|(stat)±g:S) ■ 10 6 /cm 2 s (57) 

in excellent agreement with SSM results. 

Furthermore the results point to a unique solution in the MSW plane, the 
LMA solution, with a best fit Sm 2 ~ 5 • 10~ 5 eV 2 and sin 2 26 ~ 0.75. Thus the 
oscillation corresponds to a new Sm 2 , distinct from the atmospheric v solution, 
and while the oscillation appears to be not quite maximal, the mixing angle is 
again large. For three light neutrinos, we can label this mixing as that between 
mass eigenstates 1 and 2 with 6*12 ~ 30°. 

The SuperKamiokande and SNO determinations of #23 and 6\2 are an 
exciting and important step in defining the neutrino mixing matrix, a quantity 
we hope will point the way to the next standard model. The three-flavor 
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Figure 12: The B solar neutrino flux decomposed into electron and heavy-flavor components. 
The diagonal bands show the total flux measured by the SNO neutral-current reaction (solid) 
and predicted by the SSM (dashed). The intersection with charge-current and neutral current 
bands determines the flavor content. From Ref. [8]. 



neutrino mixing matrix is conventionally written as 

C12C13 S12C13 si 3 e" l,5 \ /*t\ 

-S12C23 - Cl2S23Sl3e l<5 C12C23 - Sl2S23Sl3e" S S23C13 ^2 
S12S23 - Cl2C23Sl3e l5 -C12S23 - Sl2C23Sl3e l5 C23C13 / \ ^3 / 

1 \ / ci 3 si 3 e~ l,5 \ / C12 S12 \ f v i\ 

C23 S23 1 -sia C12 ^2 (58) 

-S23 C23 / \-si3e 115 C13 / \ 1/ \V 3 J 

Here C12 = cosf?i2, etc. Despite the discoveries of the past five years, we have 
yet to complete this matrix. Mixing between eigenstates 1 and 3 has not yet 
measured: disappearance results from the Chooz 33 and Palo Verde 34 reactor 
experiments limit #13 < 10° in the atmospheric dm 2 range. Furthermore there 
is great interest in measuring CP violation effects, due to the large mixing an- 
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Figure 13: Allowed regions of the MSW plane consistent with the SNO day and night energy 
spectra, the results from the CI and SAGE/GALLEX/GNO experiments, and the day and 
night spectra from SuperKamiokande. The star indicates the best fit. Only the LMA region 
is allowed at 99% c.l. From Ref. [32]. 



gles so far measured and to the possibility that the baryon number asymmetry 
arises through leptogenesis. CP violation in oscillations requires a nonvanish- 
ing sin #13 in addition to 5. Several very long baseline oscillation proposals have 
been made in which CP-violation would be disentangled from matter effects 
and from various neutrino parameter uncertainties. 

Recently two important laboratory neutrino oscillation results have added 
to our picture of the mixing matrix. The atmospheric neutrino 5m 2 range 
has been tested in the K2K experiment, in which events initiated by ~ 1 
GeV v^s produced by the KEK proton synchrontron are recorded in Su- 
perKamiokande 35 . The null hypothesis of no oscillations is allowed only at a 
confidence level < 0.007. The best-fit oscillation parameters, Sm 2 ~ 2.7 • 10~ 3 
eV 2 and sin 2 26 ~ 1, are in excellent agreement with the atmospheric neutrino 
values. The current data represents about half of the beam time expected for 
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Figure 14: Two-neutrino oscillation constraints from KamLAND superimposed on the solar 
neutrino LMA allowed region. From Ref. [36]. 



the experiment. 

Very recently KamLAND, the first terrestrial experiment to probe the solar 
5m 2 range, announced initial results 36 . This long-baseline experiment records 
reactor v e events in a liquid scintillator detector built at the site that once 
housed Kamiokande. Initial results confirm oscillations at 99.95% c.L, rule 
out all two-flavor solar neutrino solutions other than LMA, and significantly 
narrow the range of allowed 5m\ 2 when combined with the solar neutrino 
results. The KamLAND results are shown in fig. 14. Ultimately results from 
KamLAND combined with a high-precision measurement of solar pp neutrinos 
could determine both (5mf 2 and sin 2 29\2 with increased accuracy. 

Important question remain about the masses as well. All of the oscillation 
results test mass differences. Thus the absolute scale of v mass must be probed 
in other experiments, such as tritium /3 decay. The current limit from the 
Mainz and Troitsk experiments is m Pe < 2.2 eV 37 . As results from large-scale 
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structure surveys now underway will be affected by v masses as small as 0.3 
eV, further progress must be made in direct v mass experiments to eliminate 
a potentially significant cosmological uncertainty. 

We do not know the mass hierarchy. Both normal (the solar neutrino mass 
pair light) and inverted (this pair heavy) are allowed by the data. 

We do not know the charge conjugation properties of neutrinos. Neutrinos 
are unique among standard model fermions in allowing both Dirac and Majo- 
rana mass terms, as was noted earlier in the discussion of the seesaw mecha- 
nism. These mass terms can be distinguished because Majorana masses break 
lepton number conservation, allowing neutrinoless (3(3 decay to take place. Sev- 
eral recent next-generation (3(3 decay proposals argue that current lifetime lim- 
its can be improved by factors of <~ 100, yielding sensitivities to Majorana 
masses of 10-50 milli-eV. Many mass scenarios consistent with the solar and 
atmospheric neutrino results predict neutrinoless [3(3 decay at this level. 

Thus we are a very exciting juncture. We have found the first evidence 
for physics beyond the standard model, and we know this physics has immedi- 
ate consequences beyond nuclear and particle physics, affecting particle dark 
matter, large-scale structure, and baryogenesis. The discoveries coming from 
astrophysical neutrinos are stimulating ambitious new experiments with both 
terrestrial and astrophysical beams. A great deal of information remains hid- 
den and yet is accessible to the next generation of experiments, experiments 
that will require new neutrino beams and long-basclinc mcgadetectors, as well 
as high sensitive instruments for measuring double beta decay and probing the 
low-energy portion of the solar neutrino flux. There are many reasons to hope 
that these endeavors will provide an experimental foundation for constructing 
the next standard model of subatomic physics. 

3 Supernovae, Supernova Neutrinos, and Nucleosynthesis 

Consider a massive star, in excess of 10 solar masses, burning the hydrogen in 
its core under the conditions of hydrostatic equilibrium. When the hydrogen 
is exhausted, the core contracts until the density and temperature are reached 
where 3a — ^ 12 C can take place. The He is then burned to exhaustion. This 
pattern (fuel exhaustion, contraction, and ignition of the ashes of the previous 
burning cycle) repeats several times, leading finally to the explosive burning 
of 28 Si to Fe. For a heavy star, the evolution is rapid: the star has to work 
harder to maintain itself against its own gravity, and therefore consumes its 
fuel faster. A 25 solar mass star would go through all of these cycles in about 
7 My, with the final explosive Si burning stage taking a few days. The result 
is an "onion skin" structure of the pre-collapse star in which the star's history 
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can be read by looking at the surface inward: there are concentric shells of H, 
4 Hc, 12 C, 16 and 20 Nc, and 28 Si, with Fe at the center. 

3.1 The Explosion Mechanism 36 

The source of energy for this evolution is nuclear binding energy. A plot of the 
nuclear binding energy <5 as a function of nuclear mass shows that the minimum 
is achieved at Fe. In a scale where the 12 C mass is picked as zero: 
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Once the Si burns to produce Fe, there is no further source of nuclear energy 
adequate to support the star. So as the last remnants of nuclear burning take 
place, the core is largely supported by degeneracy pressure, with the energy 
generation rate in the core being less than the stellar luminosity The core 
density is about 2 xlO 9 g/cc and the temperature is kT ~ 0.5 MeV. 

Thus the collapse that begins with the end of Si burning is not halted 
by a new burning stage, but continues. As gravity does work on the matter, 
the collapse leads to a rapid heating and compression of the matter. As the 
nucleons in Fe are bound by about 8 MeV, sufficient heating can release as and 
a few nucleons. At the same time, the electron chemical potential is increasing. 
This makes electron capture on nuclei and any free protons favorable, 

e~+p^v e + n. (59) 

Note that the chemical equilibrium condition is 

Me + M P = + {E„} ■ (60) 

Thus the fact that neutrinos are not trapped plus the rise in the electron Fermi 
surface as the density increases, lead to increased neutronization of the matter. 
The escaping neutrinos carry off energy and lepton number. Both the electron 
capture and the nuclear excitation and disassociation take energy out of the 
electron gas, which is the star's only source of support. This means that the 
collapse is very rapid. Numerical simulations find that the iron core of the star 
(<~ 1.2-1.5 solar masses) collapses at about 0.6 of the free fall velocity. 

In the early stages of the infall the v e s readily escape. But neutrinos 
are trapped when a density of ~ 10 12 g/cm 3 is reached. At this point the 
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neutrinos begin to scatter off the matter through both charged current and 
coherent neutral current processes. The neutral current neutrino scattering 
off nuclei is particularly important, as the scattering cross section is off the 
total nuclear weak charge, which is approximately the neutron number. This 
process transfers very little energy because the mass energy of the nucleus is 
so much greater than the typical energy of the neutrinos. But momentum is 
exchanged. Thus the neutrino "random walks" out of the star. When the 
neutrino mean free path becomes sufficiently short, the "trapping time" of the 
neutrino begins to exceed the time scale for the collapse to be completed. This 
occurs at a density of about 10 12 g/cm 3 , or somewhat less than 1% of nuclear 
density. After this point, the energy released by further gravitational collapse 
and the star's remaining lepton number are trapped within the star. 

If we take a neutron star of 1.4 solar masses and a radius of 10 km, an 
estimate of its binding energy is 

GM 2 

— — ~ 2.5 x 10 53 ergs. (61) 
2rt 

Thus this is roughly the trapped energy that later will be radiated in neutrinos. 

The trapped lepton fraction Yl is a crucial parameter in the explosion 
physics: a higher trapped Yl leads to a larger homologous core, a stronger 
shock wave, and easier passage of the shock wave through the outer core, 
as will be discussed below. Most of the lepton number loss of an infalling 
mass element occurs as it passes through a narrow range of densities just 
before trapping. The reasons for this are relatively simple: on dimensional 
grounds weak rates in a plasma go as T 5 , where T is the temperature. Thus 
the electron capture rapidly turns on due to heating of the matter as it falls 
toward the trapping radius, and lepton number loss is maximal just prior to 
trapping. Inelastic neutrino reactions have an important effect on these losses, 
as the coherent trapping cross section goes as E 2 and is thus least effective 
for the lowest energy neutrinos. As these neutrinos escape, inelastic reactions 
rcpopulate the low energy states, allowing the neutrino emission to continue. 

The velocity of sound in matter rises with increasing density. The inner 
homologous core, with a mass M^c ~ 0.6 — 0.9 solar masses, is that part of the 
iron core where the sound velocity exceeds the infall velocity. In this subsonic 
central region any pressure variations that may develop in the homologous core 
during infall can smooth out before the collapse is completed. As a result, the 
homologous core collapses as a unit, retaining its density profile. That is, if 
nothing were to happen to prevent it, the homologous core would collapse to 
a point. 
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The collapse of the homologous core continues until nuclear densities are 
reached. As nuclear matter is rather incompressible (~ 200 MeV/f 3 ), the 
nuclear equation of state is effective in halting the collapse: maximum densities 
of 3-4 times nuclear are reached, e.g., perhaps 6 • 10 14 g/cm 3 . The innermost 
shell of matter reaches this supernuclear density first, rebounds, sending a 
pressure wave out through the homologous core. This wave travels faster than 
the infalling matter, as the homologous core is characterized by a sound speed 
in excess of the infall speed. Subsequent shells follow. The resulting series of 
pressure waves collect near the sonic point (the edge of the homologous core). 
As this point reaches nuclear density and comes to rest, a shock wave breaks 
out and begins its traversal of the outer core. 

Initially the shock wave may carry an order of magnitude more energy 
than is needed to eject the mantle of the star (less than 10 51 ergs). But 
as the shock wave travels through the outer iron core, it heats and melts 
the iron that crosses the shock front, at a loss of <~ 8 MeV/nucleon. The 
enhanced electron capture that occurs off the free protons left in the wake of 
the shock, coupled with the sudden reduction of the neutrino opacity of the 
matter (recall o coherent ~ N 2 ), greatly accelerates neutrino emission. This is 
another energy loss. 6 The summed losses from shock wave heating and neutrino 
emission are comparable to the initial energy carried by the shock wave. Thus 
most numerical models fail to produce a successful "prompt" hydrodynamic 
explosion: the shock stalls before it reaches the outer mantle. 

Most of the attention in the past decade focused on two explosion scenarios. 
In the prompt mechanism described above, the shock wave is sufficiently strong 
to survive the passage of the outer iron core with enough energy to blow off 
the mantle of the star. The most favorable results were achieved with smaller 
stars (less than 15 solar masses) where there is less overlying iron, and with soft 
equations of state, which produce a more compact neutron star and thus lead 
to more energy release. In part because of the lepton number loss problems 
discussed earlier, now it is widely believed that this mechanism fails for all but 
unrealistically soft nuclear equations of state. 

The delayed mechanism begins with a failed hydrodynamic explosion; after 
about 0.01 seconds the shock wave stalls at a radius of 200-300 km. It exists 
in a sort of equilibrium, gaining energy from matter falling across the shock 
front, but losing energy to the heating of that material. However, after perhaps 

6 Many numerical models predict, in conjunction with this sudden decrease in opacity, a 
strong "breakout" burst of u e s in the few milliseconds required for the shock wave to travel 
from the edge of the homologous core to the neutrinosphcrc at p ~ 10 12 g/cm 3 and r ~ 50 
km. The neutrinosphcrc is the term for the neutrino trapping radius, or surface of last 
scattering. 
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0.5 seconds, the shock wave is revived due to neutrino heating of the nucleon 
"soup" left in the wake of the shock. This heating comes primarily from 
charged current reactions off the nucleons in that nucleon gas; quasi-elastic 
scattering also contributes. This high entropy radiation-dominated gas may 
reach two MeV in temperature. The pressure exerted by this gas helps to push 
the shock outward. It is important to note that there are limits to how effective 
this neutrino energy transfer can be: if matter is too far from the core, the 
coupling to neutrinos is too weak to deposit significant energy. If too close, the 
matter may be at a temperature (or soon reach a temperature) where neutrino 
emission cools the matter as fast or faster than neutrino absorption heats it. 
The term "gain radius" is used to describe the region where useful heating is 
done. 

This subject is still controversial and unclear. The problem is numerically 
challenging, forcing modelers to handle the difficult hydrodynamics of a shock 
wave; the complications of the nuclear equation of state at densities not yet 
accessible to experiment; modeling in two or three dimensions; handling the 
slow diffusion of neutrinos; etc. Not all of these aspects can be handled rea- 
sonably at the same time, even with existing supercomputers. Thus there is 
considerable disagreement about whether we have any supernova model that 
succeeds in ejecting the mantle. 

However the explosion proceeds, there is agreement that 99% of the 3 -10 53 
ergs released in the collapse is radiated in neutrinos of all flavors. The time 
scale over which the trapped neutrinos leak out of the protoncutron star is 
about three seconds. Through most of their migration out of the protoneutron 
star, the neutrinos are in flavor equilibrium 

e.g., v e + i> e <-+ i/p + Dfj,. (62) 

As a result, there is an approximate equipartition of energy among the neutrino 
flavors. After weak decoupling, the v e s and i7 e s remain in equilibrium with the 
matter for a longer period than their heavy-flavor counterparts, due to the 
larger cross sections for scattering off electrons and because of the charge- 
current reactions 

v e + n <-> p + e~ 

i7 e + p ^n + e + . (63) 

Thus the heavy flavor neutrinos decouple from deeper within the star, where 
temperatures are higher. Typical calculations yield 

T V)x ~ T Vt ~ 8MeV T„ o - 3.5MeV TV. - 4.5MeV. (64) 
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The difference between the v e and v e temperatures is a result of the neutron 
richness of the matter, which enhances the rate for charge-current reactions of 
the v e s, thereby keeping them coupled to the matter somewhat longer. (This 
temperature hierarchy, particular the difference between the heavy-flavor and 
electron neutrino temperatures, is still a matter of some debate, as it is in- 
fluenced by both the explosion mechanism (e.g., by neutron fingers and other 
types of mixing) and the detailed modeling of the microphysics. Some recent 
work 39 argues for considerably smaller temperature differences than those 
given above.) 

This temperature hierarchy is important because temperature inversions 
are a potential signature of oscillations. A three-flavor MSW level-crossing 
diagram is shown in fig. 15. Naively (that is, without considering neutrino- 
neutrino scattering and other effects that could alter this picture) the crossings 
corresponding to 5ra 2 solar and Sm^ tmos both occur well outside the neutrino 
sphere, that is, after the neutrinos have decoupled and have fixed spectra 
characterized by the temperatures given above. (For example, 5m 2 atmos corre- 
sponds to a electron density typical of the base of the carbon zone, prior to the 
explosion.) Thus a v e <-» v T oscillation would produce a distinctive T ~ 8 MeV 
spectrum of v e s. Because neutrino-nucleus cross sections often grow as a high 
power of the neutrino energy (due both to phase space and threshold effects), 
this will produce elevated v e event rates in many detectors. Oscillations may 
also have an effect on nucleosynthesis, such as the ^-process we will discuss 
below. 

3.2 The Neutrino Process 40 

Core-collapse supernovae are one of the major engines driving galactic chem- 
ical evolution, producing and ejecting the metals that enrich our galaxy. The 
discussion of the previous section described the hydrostatic evolution of a pre- 
supernova star in which large quantities of the most abundant metals (C, O, 
Ne,...) are synthesized and later ejected during the explosion. During the pas- 
sage of the shock wave through the star's mantle, temperature of ~ (1— 3)T0 9 K 
and are reached in the silicon, oxygen, and neon shells. This shock wave heat- 
ing induces (7, a) <-> (a, 7) and related reactions that generate a mass flow 
toward highly bound nuclei, resulting in the synthesis of iron peak elements 
as well as less abundant odd-A species. Rapid neutron-induced reactions are 
thought to take place in the high-entropy atmosphere just above the mass cut, 
producing about half of the heavy elements above A <~ 80. Finally, the v- 
process described below is responsible for the synthesis of rare species such as 
n B and 19 F. This process involves the weak response of nuclei at momentum 
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~10 g/cm density vacuum 

Figure 15: Three-flavor neutrino level-crossing diagram. The illustrated u e <-> and v e <-> 
v T crossings should arise from the solar and atmospheric Sm 2 s. The diagram illustrates 
that these crossing occur outside the neutrinosphcrc, in regions where the neutrinos are fully 
decoupled from the matter. 

transfers where the allowed approximation is no longer valid. Thus we will use 
the ^-process in this section to illustrate some of the relevant nuclear physics. 

One of the problems - still controversial - that may be connected with the 
neutrino process is the origin of the light elements Be, B and Li, elements which 
are not produced in sufficient amounts in the big bang or in any of the stellar 
mechanisms we have discussed. The traditional explanation has been cosmic 
ray spallation interactions with C, O, and N in the interstellar medium. In this 
picture, cosmic ray protons collide with C at relatively high energy, knocking 
the nucleus apart. So in the debris one can find nuclei like 10 B, n B, and 7 Li. 

But there are some problems with this picture. First of all, this is an 
example of a secondary mechanism: the interstellar medium must be enriched 
in the C, O, and N to provide the targets for these reactions. Thus cosmic ray 
spallation must become more effective as the galaxy ages. The abundance of 
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boron, for example, would tend to grow quadratically with metallicity, since 
the rate of production goes linearly with metallicity. But observations, espe- 
cially recent measurements with the HST, find a linear growth 41 in the boron 
abundance. 

A second problem is that the spectrum of cosmic ray protons peaks near 
1 GeV, leading to roughly comparable production of the two isotopes 10 B 
and 11 B. That is, while it takes more energy to knock two nucleons out of 
carbon than one, this difference is not significant compared to typical cosmic 
ray energies. More careful studies lead to the expectation that the abundance 
ratio of n B to 10 B might be <~ 2. In nature, it is greater than 4. 

Fans of cosmic ray spallation have offered solutions to these problems, 
e.g., similar reactions occurring in the atmospheres of nebulae involving lower 
energy cosmic rays. As this suggestion was originally stimulated by the ob- 
servation of nuclear 7 rays from Orion, now retracted, some of the motivation 
for this scenario has evaporated. Here we focus on an alternative explanation, 
synthesis via neutrino spallation. 

Previously we described the allowed Gamow- Teller (spin-flip) and Fermi 
weak interaction operators. These are the appropriate operators when one 
probes the nucleus at a wavelength that is, at a size scale - where the nucleus 
responds like an elementary particle. We can then characterize its response by 
its macroscopic quantum numbers, the spin and charge. On the other hand, 
the nucleus is a composite object and, therefore, if it is probed at shorter length 
scales, all kinds of interesting radial excitations will result, analogous to the 
vibrations of a drumhead. For a reaction like neutrino scattering off a nucleus, 
the full operator involves the additional factor 

e lt? ~ 1 + ik ■ r (65) 

where the expression on the right is valid if the magnitude of k is not too large. 
Thus the full charge operator includes a "first forbidden" term 

A 

E r >3« (66) 

i=i 

and similarly for the spin operator 

A 

$^Fi®*(*)]j=0,l,27S(*)- (67) 
i=l 

These operators generate collective radial excitations, leading to the so-called 
"giant resonance" excitations in nuclei. The giant resonances are typically at 
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an excitation energy of 20-25 MeV in light nuclei. One important property is 
that these operators satisfy a sum rule (Thomas-Reiche-Kuhn) of the form 



where the sum extends over a complete set of final nuclear states. These 
first-forbidden operators tend to dominate the cross sections for scattering the 
high energy supernova neutrinos (u^s and v r s), with E v <~ 25 MeV, off light 
nuclei. From the sum rule above, it follows that nuclear cross sections per 
target nucleon are roughly constant. 

The El giant dipole mode described above is depicted qualitatively in fig. 
16a. This description, which corresponds to an early model of the giant res- 
onance response by Goldhaber and Teller, involves the harmonic oscillation 
of the proton and neutron fluids against one another. The restoring force for 
small displacements would be linear in the displacement and dependent on the 
nuclear symmetry energy. There is a natural extension of this model to weak 
interactions, where axial excitations occur. For example, one can envision a 
mode similar to that of fig. 16a where the spin-up neutrons and spin-down pro- 
tons oscillate against spin-down neutrons and spin-up protons, the spin-isospin 
mode of fig. 16b. This mode is one that arises in a simple SU(4) extension of 
the Goldhaber- Teller model, derived by assuming that the nuclear force is spin 
and isospin independent, at the same excitation energy as the El mode. In 
full, the Goldhaber- Teller model predicts a degenerate 15-dimcnsional super- 
multiplet of giant resonances, each obeying sum rules analogous to the TRK 
sum rule. While more sophisticated descriptions of the giant resonance region 
are available, of course, this crude picture is qualitatively accurate. 

This nuclear physics is important to the ^-process. The simplest example 
of ^-process nucleosynthesis involves the Ne shell in a supernova. Because of 
the first-forbidden contributions, the cross section for inelastic neutrino scat- 
tering to the giant resonances in Ne is <~ 3 • 10~ 41 cm 2 /fiavor for the more 
energetic heavy-flavor neutrinos. This reaction 



transfers an energy typical of giant resonances, ~ 20 MeV. A supernova releases 
about 3 xlO 53 ergs in neutrinos, which converts to about 4 x 10 57 heavy flavor 
neutrinos. The Ne shell in a 20 M Q star has at a radius ~ 20,000 km. Thus 
the neutrino fluence through the Ne shell is 




(68) 



v + A -y v' + A* 



(69) 



4 • 10 57 



10 38 /cm 2 . 
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Figure 16: Schematic illustration of a) the El giant dipole mode familiar from electromag- 
netic interactions and b) a spin-isospin giant dipole mode associated with the first-forbidden 
weak axial response. 



Thus folding the fluence and cross section, one concludes that approximately 
l/300th of the Ne nuclei interact. 

This is quite interesting since the astrophysical origin of 19 F had not been 
understood. The only stable isotope of fluorine, 19 F has an abundance 



19t 



20 Ne 3100' 



(71) 



This leads to the conclusion that the fluorine found in a tube of toothpaste 
was created by neutral current neutrino reactions deep inside some ancient 
supernova. 

The calculation of the final 19 F/ 20 Ne ratio is more complicated than the 
simple 1/300 ratio given above: 

• When Ne is excited by ~ 20 MeV through inelastic neutrino scattering, it 
breaks up in two ways 
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Ne(i/, i/) 20 Ne* 
Ne(i/, z/) 20 Ne* 



> 19 Nc + n 
> 19 F + p 
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F + c+ + v c + n 



(72) 



with the first reaction occurring half as frequently as the second. As both 
channels lead to 19 F, we have correctly estimated the instantaneous abundance 
ratio in the Ne shell of 

19 F 1 

2^ ~ 300" (?3) 



4G 



• We must also address the issue of whether the produced 19 F survives. In the 
first 10 -8 sec the co-produced neutrons in the first reaction react via 

15 0(n,p) 15 N 19 Nc(n,a) 16 20 Nc(n, 7 ) 21 Ne 19 Ne(n,p) 19 F (74) 

with the result that about 70% of the 19 F produced via spallation of neutrons 
is then immediate destroyed, primarily by the (n, a) reaction above. In the 
next 10 -6 sec the co-produced protons are also processed 

15 N(p,a) 12 C 19 F(p,a) 16 23 Na(p, a) 20 Ne (75) 

with the latter two reactions competing as the primary proton poisons. This 
makes an important prediction: stars with high Na abundances should make 
more F, as the 23 Na acts as a proton poison to preserve the produced F. • 
Finally, there is one other destruction mechanism, the heating associated with 
the passage of the shock wave. It turns out the the F produced prior to shock 
wave passage can survive if it is in the outside half of the Ne shell. The reaction 

19 F( 7 ,a) 15 N (76) 

destroys F for peak explosion temperatures exceeding 1.7-10 9 K. Such a temper- 
ature is produced at the inner edge of the Ne shell by the shock wave heating, 
but not at the outer edge. 

If all of this physics in handled is a careful network code that includes 
the shock wave heating and F production both before and after shock wave 
passage, the following are the results: 

[ 19 F/ 20 Ne]/[ 19 F/ 20 Ne] Q T hcavy „(MeV) 



0.14 4 

0.6 6 

1.2 8 

1.1 10 

1.1 12 



where the abundance ratio in the first column has been normalized to the solar 
value. One sees that the attribution of F to the neutrino process argues that 
the heavy flavor v temperature must be greater than 6 MeV, a result theory 
favors. One also sees that F cannot be overproduced by this mechanism: 
although the instantaneous production of F continues to grow rapidly with the 
neutrino temperature, too much F results in its destruction through the (p, a) 
reaction, given the metalicity assumed in this calculation (a solar abundance 
of the competing proton poison 23 Na). Indeed, this illustrates an odd quirk: 
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although in most cases the neutrino process is a primary mechanism, one needs 
23 Na present to produce significant F. Thus in this case the neutrino process 
is a secondary mechanism. 

While there are other significant neutrino process products ( 7 Li, 138 La, 
180 Ta, 15 N ...), the most important product is n B, produced by spallation off 
carbon. A calculation by Timmes et al. 41 found that the combination of the 
neutrino process, cosmic ray spallation and big-bang nucleosynthesis together 
can explain the evolution of the light elements. The neutrino process, which 
produces a great deal of 11 B but relatively little 10 B, combines with the cos- 
mic ray spallation mechanism to yield the observed isotope ratio. Again, one 
prediction of this picture is that early stars should be 11 B rich, as the neutrino 
process is primary and operates early in our galaxy's history; the cosmic ray 
production of 10 B is more recent. There is hope that HST studies will soon be 
able to discriminate between 10 B and 11 B: as yet this has not been done. 



3.3 The r-process 

Beyond the iron peak nuclear Coulomb barriers become so high that charged 
particle reactions become ineffective, leaving neutron capture as the mechanism 
responsible for producing the heaviest nuclei. If the neutron abundance is 
modest, this capture occurs in such a way that each newly synthesized nucleus 
has the opportunity to (3 decay, if it is energetically favorable to do so. Thus 
weak equilibrium is maintained within the nucleus, so that synthesis is along 
the path of stable nuclei. This is called the s- or slow-process. However a plot 
of the s-process in the (N,Z) plane reveals that this path misses many stable, 
neutron-rich nuclei that are known to exist in nature. This suggests that 
another mechanism is at work, too. Furthermore, the abundance peaks found 
in nature near masses A ~ 130 and A ~ 190, which mark the closed neutron 
shells where neutron capture rates and (3 decay rates are slower, each split into 
two sub-peaks. One set of sub-peaks corresponds to the closed-neutron-shell 
numbers N <~ 82 and N ~ 126, and is clearly associated with the s-process. The 
other set is shifted to smaller N, <~ 76 and <~ 116, respectively, and is suggestive 
of a much more explosive neutron capture environment where neutron capture 
can be rapid. 

This second process is the r- or rapid-process, characterized by: 

• The neutron capture is fast compared to (3 decay rates. 

• The equilibrium maintained within a nucleus is established by (n, 7) ^> (7,n): 
neutron capture fills up the available bound levels in the nucleus until this 
equilibrium sets in. The new Fermi level depends on the temperature and the 
relative n/7 abundance. 
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• The nucleosynthesis rate is thus controlled by the decay rate: each 0~ 
capture converting n — > p opens up a hole in the neutron Fermi sea, allowing 
another neutron to be captured. 

• The nucleosynthesis path is along exotic, neutron-rich nuclei that would be 
highly unstable under normal laboratory conditions. 

• As the nucleosynthesis rate is controlled by the decay, mass will build up 
at nuclei where the decay rates are slow. It follows, if the neutron flux is 
reasonable steady over time so that equilibrated mass flow is reached, that the 
resulting abundances should be inversely proportional to these decay rates. 

Let's first explore the (n, 7) <-» (7,11) equilibrium condition, which requires 
that the rate for (n, 7) balances that for (7,11) for an average nucleus. So 
consider the formation cross section 



A + n^(A + l)+ 7 . (77) 

This is an exothermic reaction, as the neutron drops into the nuclear well. Our 
averaged cross section, assuming a resonant reaction (the level density is high 
in heavy nuclei) is 

<"><«» = {wY T ^ ,KT (78) 

where E ~ is the resonance energy, and the Ts are the indicated partial and 
total widths. Thus the rate per unit volume is 

where N n and Na are the neutron and nuclear number densities and \x the 
reduced mass. This has to be compared to the (7,11) rate. 

The (7,n) reaction requires the photon number density in the gas. This is 
given by the Bose-Einstein distribution 

at/ \ 87r e 2 de /0 „. 

N ( € ) = 575 iwr 7- 80 

c 3 ft, 3 e £ / fcT - 1 

The high-energy tail of the normalized distribution can thus be written 



< 81 > 

where in the last expression we have set % = c = 1. 
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Now we need the resonant cross section in the (7,n) direction. For photons 
the wave number is proportional to the energy, so 



As the velocity is c =1, 

M ^f eV ' ,1T *5(^T(W' (83) 

We evaluate this in the usual way for a sharp resonance, remembering that the 
energy integral over just the denominator above (the sharply varying part) is 
2tt/T 

r r ......■> 

(84) 



r 7 r n E r /kT 2 



7V 7 



So that the rate becomes 



r (7in) ~2AT A+1 ^e-^/ feT . (85) 
Equating the (n, 7) and (7,n) rates and taking Na ~ Na-i then yields 



{he) 3 V 2tt 

where the hs and cs have been properly inserted to give the right dimensions. 
Now E r is essentially the binding energy. So plugging in the conditions N n ~ 
3 x 10 23 /cm 3 and T 9 ~ 1, we find that the binding energy is ~ 2.4 MeV. Thus 
neutrons are bound by about 30 times kT, a value that is still small compared 
to a typical binding of 8 MeV for a normal nucleus. (In this calculation the 
neutron reduced mass is calculated by assuming a nuclear target with A=150.) 

The above calculation fails to count spin states for the photons and nuclei 
and is thus not quite correct. But it makes the essential point: the r-process 
involves very exotic species largely unstudied in any terrestrial laboratory. It 
is good to bear this in mind, as in the following section we will discuss the 
responses of such nuclei to neutrinos. Such responses thus depend on the 
ability of theory to extrapolate responses from known nuclei to those quite 
unfamiliar. 

The path of the r-process is along neutron-rich nuclei, where the neutron 
Fermi sea is just ~ (2-3) MeV away from the neutron drip line (where no 
more bound neutron levels exist). After the r-process finishes (the neutron 
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exposure ends) the nuclei decay back to the valley of stability by /3 decay. This 
can involve some neutron spallation (/3-delayed neutrons) that shift the mass 
number A to a lower value. But it certainly involves conversion of neutrons 
into protons, and that shifts the r-process peaks at N ~ 82 and 126 to a lower 
N, off course. This effect is clearly seen in the abundance distribution: the 
r-process peaks are shifted to lower N relative to the s-process peaks. This 
is the origin of the second set of "sub-peaks" mentioned at the start of the 
section. 

It is believed that the r-process can proceed to very heavy nuclei (A ~ 
270) where it is finally ended by /3-delayed and n-induced fission, which feeds 
matter back into the process at an A <~ A max /2. Thus there may be important 
cycling effects in the upper half of the r-process distribution. 

What is the site(s) of the r-process? This has been debated many years 
and still remains a controversial subject: 

• The r-process requires exceptionally explosive conditions 

p(n) - 10 20 cm" 3 T - 10 9 K t - Is. 

• Both have been primary and secondary sites have been proposed. Primary 
sites are those not requiring preexisting metals. Secondary sites are those 
where the neutron capture occurs on preexisting s-process seeds. 

• Suggested primary sites include the the neutronized atmosphere above the 
proto-neutron star in a Type II supernova, neutron-rich jets produced in su- 
pernova explosions or in neutron star mergers, inhomogeneous big bangs, etc. 

• Secondary sites, where p(n) can be lower for successful synthesis, include the 
He and C zones in Type II supernovae, the red giant He flash, etc. 

The balance of evidence favors a primary site, so one requiring no pre- 
enrichment of heavy s-process metals. Among the evidence: 

1) Keck and HST studies of very- metal-poor halo stars: The most important 
evidence are the recent measurements of Cowan, Sneden et al. 43 of very metal- 
poor stars ([Fe/H] ~ -1.7 to -3.12) where an r-process distribution very much 
like that of our sun has been seen for Z > 56. Furthermore, in these stars 
the iron content is variable. This suggests that the "time resolution" inherent 
in these old stars is short compared to galactic mixing times (otherwise Fe 
would be more constant). The conclusion is that the r-process material in 
these stars is most likely from one or a few local supernovae. The fact that 
the distributions match the solar r-process (at least above charge 56) strongly 
suggests that there is some kind of unique site for the high-Z portion of the 
r-process: the solar r-process distribution did not come from averaging over 
many different kinds of r-process events. Clearly the fact that these old stars 
are enriched in r-process metals also strongly argues for a primary process: 
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the r-process works quite well in an environment where there are few initial s- 
process metals. (The situation for elements below Z=56 is somewhat different. 
While an adequate discussion would take us beyond the limits of these lectures, 
interested reader are directed to Ref. [42].) 

2) There are also fairly good theoretical arguments that a primary r-process oc- 
curring in a core-collapse supernova might be viable . First, galactic chemical 
evolution studies indicate that the growth of r-process elements in the galaxy 
is consistent with low-mass Type II supernovae in rate and distribution. More 
convincing is the fact that modelers have shown that the conditions needed for 
an r-process (very high neutron densities, temperatures of 1-3 billion degrees) 
might be realized in a supernova. The site is the last material expelled from 
the supernova, the matter just above the mass cut. When this material is 
blown off the star initially, it is a very hot neutron-rich, radiation-dominated 
gas containing neutrons and protons, but an excess of the neutrons. As it 
expands off the star and cools, the material first goes through a freeze-out to 
a particles, a step that essentially locks up all the protons in this way. Then 
the as interact through reactions like 

a + a + a —> 12 C 
a + a + n — > 9 Be 

to start forming heavier nuclei. Note, unlike the big bang, that the density 
is high enough to allow such three-body interactions to bridge the mass gaps 
at A = 5,8. The a capture continues up to heavy nuclei, to A ~ 80, in the 
network calculations. The result is a small number of "seed" nuclei, a large 
number of as, and excess neutrons. These neutrons preferentially capture on 
the heavy seeds to produce an r-process. Of course, what is necessary is to 
have ~ 100 excess neutrons per seed in order to successfully synthesize heavy 
mass nuclei. Some of the modelers find conditions where this almost happens. 

There are some very nice aspects of this site: the amount of matter ejected 
is about 10~ 5 — 10~ 6 solar masses, which is just about what is needed over the 
lifetime of the galaxy to give the integrated r-process metals we see, taking a 
reasonable supernova rate. But there are also a few problems, especially the 
fact that with calculated entropies in the nucleon soup above the proto-neutron 
star, neutron fractions appear to be too low to produce a successful A <~ 190 
peak. There is some interesting recent work invoking neutrino oscillations 45 
to cure this problem: charge current reactions on free protons and neutrons 
determine the n/p ratio in the gas. This is discussed briefly in the next section. 

The nuclear physics of the r-process tells us that the synthesis occurs 
when the nucleon soup is in the temperature range of (3-1) T0 9 K, which, in 
the hot bubble r-process described above, corresponds to a freeze-out radius 
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of (600-1000) km and a time ~ 10 seconds after core collapse. The neutrino 
fluence after freeze-out (when the temperature has dropped below 10 9 K and 
the r-process stops) is then ~ (0.045-0.015) -10 51 ergs/(100km). Thus, after 
completion of the r-process, the newly synthesized material experiences an 
intense flux of neutrinos. This brings up the question of whether the neutrino 
flux could have any effect on the r-process. 

3.4 Neutrinos and the r-process 46 

Before describing the exotic effects of neutrino oscillations on the supernovae 
and nucleosynthesis, we will examine standard-model effects that are neverthe- 
less quite interesting. The nuclear physics of this section - neutrino-induced 
neutron spallation reactions - is also relevant to recently proposed supernova 
neutrino observatories such as OMNIS and LAND. Comparing to our earlier 
discussion of carbon- and neon-zone synthesis by the ^-process, it is apparent 
that neutrino effects could be much larger in the hot bubble r-process: the 
synthesis occurs much closer to the star than our Ne radius of 20,000 km: 
estimates are 600-1000 km. The r-process is completed in about 10 seconds 
(when the temperature drops to about one billion degrees), but the neutrino 
flux is still significant as the r-process freezes out. The net result is that the 
"post-processing" neutrino fluence - the fluence that can alter the nuclear dis- 
tribution after the r-process is completed - is about 100 times larger than that 
responsible for fluorine production in the Ne zone. Recalling that 1/300 of 
the nuclei in the Ne zone interacted with neutrinos, and remembering that the 
relevant neutrino- nucleus cross sections scale as A, one quickly sees that the 
probability of a r-process nucleus interacting with the neutrino flux is approx- 
imately unity. 

Because the hydrodynamic conditions of the r-process are highly uncer- 
tain, one way to attack this problem is to work backward in time. We know 
the final r-process distribution (what nature gives us) and we can calculate 
neutrino-nucleus interactions relatively well. Thus from the observed r-process 
distribution (including neutrino post-processing) we can work backward to find 
out what the r-process distribution looked like at the point of freeze-out. In 
Figs. 17 and 18, the "real" r-process distribution - that produced at freeze-out 
- is given by the dashed lines, while the solid lines show the effects of the neu- 
trino post-processing for a particular choice of fluence. The nuclear physics 
input into these calculations is precisely that previously described: GT and 
first-forbidden cross sections, with the responses centered at excitation ener- 
gies consistent with those found in ordinary, stable nuclei, taking into account 
the observed dependence on \N — Z\. 
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Figure 17: Comparison of the r-process distribution that would result from the freeze-out 
abundances near the A ~ 130 mass peak (dashed line) to that where the effects of neutrino 
post-processing have been include (solid line). The fiuence has been fixed by assuming that 
the A = 124-126 abundances are entirely due to the ^-process. 

One important aspect of the figures is that the mass shift is significant. 
This has to do with the fact that a 20 MeV excitation of a neutron-rich nucleus 
allows multiple neutrons ( <~ 5) to be emitted. (Remember we found that the 
binding energy of the last neutron in an r-process neutron-rich nuclei was about 
2-3 MeV under typical r-process conditions.) The second thing to notice is that 
the relative contribution of the neutrino process is particularly important in 
the "valleys" beneath the mass peaks: the reason is that the parents on the 
mass peak are abundant, and the valley daughters rare. In fact, it follows 
from this that the neutrino process effects can be dominant for precisely seven 
isotopes (Te, Re, etc.) lying in these valleys. Furthermore if an appropriate 
neutrino fluence is picked, these isotope abundances are produced perfectly 
(given the abundance errors). The fluences are 

N = 82 peak 0.031 • 10 51 ergs/(100km) 2 /navor 
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N = 126 peak 0.015 • 10 51 crgs/(100km) 2 /fiavor, 

values in fine agreement with those that would be found in a hot bubble r- 
process. So this is circumstantial but significant evidence that the material 
near the mass cut of a Type II supernova is the site of the r-process: there is 
a neutrino fingerprint. 




170 175 180 185 190 195 200 
A 

Figure 18: As in fig. 17, but for the A ~ 195 mass peak. The A = 183-187 abundances are 
entirely attributed to the ^-process. 

4 Neutrino Oscillation Constraints from the r-process 

Many aspects of supernova physics could be altered by neutrino oscillations. 
One satisfying consequence of the recent solar and atmospheric neutrino os- 
cillation discoveries is that the derived parameters, the mass differences and 
mixing angles, have begun to narrow some of "parameter space" for such su- 
pernova oscillation effects. 

Earlier it was noted that understanding neutrino transport is an essen- 
tial but highly nontrivial part of the supernova mechanism. The addition of 
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neutrino oscillations adds another layer of complication, and thus opportunity 
for surprises. To set the stage for this discussion, we recall that the density 
at the neutrinosphere is ~ 10 12 g cm~ 3 and the density at the position of the 
stalled shock is 47 ~2x 10 7 g cm~ 3 . Writing the MSW resonance density in 
appropriate units 

o, ,«« [5m 2 \ /10MeV\ /0.5\ nn , 
Prcs = 1.31 x 10 6 ^— — j i^—j cos 29 gem- 3 , (87) 

one sees, for E v <~ 10 MeV, that the range of 5m 2 s producing crossings between 
the neutrinosphere and the stalled shock wave is ~ 1 — 10 4 cV 2 . Some time 
ago Fuller et al. 47 pointed out that v e — v heavy oscillations in this region could 
substantial increase the explosion energy, as hot v e s would couple more strongly 
to the matter because of the increase in the v e + n cross section (which goes as 
E 2 , roughly) . Unfortunately the necessary range of Sm 2 includes neither the 
atmospheric nor solar values. Thus at least in the simplest scenario of three 
light neutrinos, the oscillations we know about from experiment may not affect 
this aspect of the explosion mechanism. (However there are many scenarios 
- a fourth neutrino motivated by LSND 48 , CPT violation, etc - where this 
conclusion would need to be reexamined. Even in the simplest three-flavor 
scenario discussed above, one should consider complications that might move 
the MSW resonance position, such as the contributions of v — v scattering to 
the MSW potential.) 

The potential for neutrino oscillation to alter the r-process (if core-collapse 
supernovae are indeed the site) was also recognized some time ago. The re- 
process requires a neutron-rich environment: the ratio of electrons to baryons, 
Y" e , should be less than one half. Y e in the nucleosynthesis region is given 
approximately 49 by 

1 + Aj7 eP / A„ e „ 1 + Tn e /T Ue ' 

where A„ e „, etc. are the capture rates. Hence if T Ve > T Ve , then the medium 
is neutron rich. But oscillations of the type v e <-» Vheavy could invert this 
hierarchy, thus destroying a necessary condition for the r-process. In the sim- 
plest scenario of three light neutrinos and 5m 2 consistent with the solar and 
atmospheric v results, this catastrophe is avoided for the hot-bubble r-process 
because the crossings occur outside the region of interest. (However, the same 
caveats noted above apply.) 

If the supernova is an r-process site it is also desirable to have a neutron to 
seed-nucleus ratio > 100 in order that the heavier r-process species (i.e., those 
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in the A = 195 peak) can be produced. This ratio is basically determined by 
three quantities: i) the expansion rate; ii) the electron fraction Y e ; and iii) the 
entropy per baryon. Though different calculations 44 ' 50 disagree on the value of 
the entropy in the neutrino-driven wind during the r-process nucleosynthesis, 
several models can produce values of these three parameters adequate for a 
successful r-process - that is, until certain v reactions are considered. These 
v reactions, acting during or immediately after freezeout, allow too many seed 
nuclei to form, thus producing an unfavorable neutron/seed ratio. This pre- 
vents the r-process from synthesizing the heaviest nuclei, those in the A <~ 190 
peak. The worst of these v reactions is called the alpha effect. 

The alpha effect 51 occurs at the epoch of alpha particle formation. As the 
temperature drops, essentially all of the protons and most of the neutrons in 
the ejecta lock themselves into tightly bound alpha particles. As the matter 
was initially neutron rich, free neutrons remain and Y e is below 0.5. The alpha 
effect pushes the electron fraction higher, towards Y e = 0.5. The increase in 
Y e comes about because v e s capture on neutrons, with the produced protons 
capturing more neutrons to produce as, thereby reducing the number of free 
neutrons available for the r-process 52 . This effect has been shown to be the 
biggest impediment to achieving an acceptable r-process yield 53 . Matter- 
enhanced neutrino oscillations between electron neutrinos and other active 
species worsen this problem, tending to increase electron neutrino energies and 
thus the average charge current cross section. On the other hand, active-sterile 
oscillations which diminish the flux of electron neutrinos outside some radius 
can reduce the alpha effect, and thus preserve the r-process. This scenario has 
been discussed extensively because sterile neutrinos may also account for the 
LSND results. 

There is an important constraint on this "solution," however. In models 
where the r-process material is blown off the protoneutron star by a neutrino 
wind, a large flux of electron neutrinos is essential to overcome the binding 
effects of gravity. As nucleons are gravitationally bound by about ~ 100 MeV 
near the surface of the protoneutron star while each neutrino has an energy 
<~ 10 MeV, <~ 10 charge-current i/-nucleon interactions arc needed to eject 
the nucleon to infinity. This in turn requires that the active-sterile oscillation 
occur at a relatively large radius, so that strong wind effects at smaller r 
are unaffected. It proves possible to arrange the necessary condition through 
active-sterile v e ^ v s and v e ^ v s channels. In such a scheme 45 the lightest 
sterile neutrino would be heavier than the v e and split from it by a vacuum 
mass-squared difference of 3 eV 2 < 5m 2 es < 70 cV 2 with vacuum mixing angle 
sin 2 28 es > 10~ 4 . Whether this solution is necessary, though, is quite another 
matter. It has been argued that very fast expansion rates could circumvent 
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the alpha effect by minimizing the v fluence experienced by the matter. This 
underscores how difficult it is to assess our current understanding of the r- 
process when our theoretical models of supernovae seem to be less successful 
than Nature in producing explosions. (This reminds one that the foundation 
for our discoveries in solar neutrino physics was a reliable SSM!) 

Finally, regardless of where the crossings occur within a supernova, the 
supernovae neutrinos that SNO, SuperKamiokande, and other detectors record 
will have been altered by oscillations, including potentially oscillation channels 
we have not yet probed elsewhere (e.g., #13). The extent to which we can 
exploit neutrinos as a probe of supernova physics will depend on how well we 
succeed with current efforts to determine the entries in the mixing matrix. 

This work was supported in part by the US Department of Energy under 
grants DE-FG03-00ER41132 and DE-FG02-01ER41187. Portions of these lec- 
tures were based on summer school notes prepared in collaboration with A. 
Baha Balantekin, whom I thank. 

References 

1. This section is an updated version of W. C. Haxton, Ann. Rev. Astron. 
Astrophys. 33, 459 (1995); the historical discussion is given in more 
detail in W. C. Haxton and B. R. Holstcin, Am. J. Phys. 68, 15 (2000). 

2. R. Davis, Jr., D. S. Harmer, and K. C. Hoffman, Phys. Rev. Lett. 20, 
1205 (1966). 

3. J. N. Abdurashitov et al., Phys. Lett. B 328, 234 (1994) and JETP 95, 
181 (2002). 

4. P. Ansclmann ct al., Phys. Lett. B 285, 376 (1992). 

5. GNO Collaboration, M. Altmann et al., Phys. Lett. B 490, 16 (2002). 

6. Y. Suzuki, Nucl. Phys. B 38, 54 (1995). 

7. Y. Suzuki, talk presented at Neutrino '98 (Takayama, Japan, June, 1998). 

8. SNO Collaboration, Phys. Rev. Lett. 89, 011301 (2002) and Phys. Rev. 
Lett. 87, 071301 (2001). 

9. Y. Fukuda et al. (Supcrkamiokande collaboration), Phys. Rev. Lett. 81, 
1562 (1998); Phys. Lett. B 433, 9 (1998). 

10. J. N. Bahcall, S. Basu, and M. H. Pinsonneault, Phys. Lett. B 433, 1 
(1998) and Ap. J. 555, 990 (2001). 

11. A. S. Brun, S. Turck-Chieze, and P. Morel, Ap. J. 506, 913 (1998) and 
private communication; S. Turck-Chieze and I. Lopez, Ap. J. 408, 347 
(1993). 

12. This section taken from W. C. Haxton, nucl-th/9901037. 

13. K. Lande, talk presented at Neutrino '98 (Takayama, Japan, June, 1998). 



58 



14. J. N. Bahcall, Neutrino Astrophysics, (Cambridge University, Cam- 
bridge, 1989). 

15. S..Fukuda et al. (Superkamiokande Collaboration), Phys. Rev. Lett. 86, 
5651 (2001); M. B. Smy, hcp-ex/0208004/ 

16. J. N. Bahcall and R. Ulrich, Rev. Mod. Phys. 60, 297 (1988). 

17. J. N. Bahcall and W. C. Haxton, Phys. Rev. D 40, 931 (1989). 

18. V. Castcllani, S. DegPInnocenti, G. Fiorentini, M. Lissia, and B. Ricci, 
Phys. Rev. D 50, 4749 (1994). 

19. K. M. Hccgcr and R. G. H. Robertson, Phys. Rev. Lett. 77, 3720 (1996). 

20. J. Christensen-Dalsgaard, Lecture Notes on Stellar Oscillations, 
http:/ / www. obs. aau. dk/~ jcd/ oscilnotes/. 

21. M. Gcll-Mann, P. Ramond, and R. Slansky, in Supergravity, P. van 
Nieuwenhuizen and D.Z. Freedman, Eds., North Holland, Amsterdam, 
1979, p. 315; T. Yanagida, in Proceedings of the Workshop on the Uni- 
fied Theory and Baryon Number in the Universe, O. Sawada and A. 
Sugamoto, Eds., (KEK Report 79-18, 1979), p.95. 

22. S. P. Mikheyev and A. Smirnov, Sov. J. Nucl. Phys. 42, 913 (1985); L. 
Wolfenstein, Phys. Rev. D 17, 2369 (1979). 

23. H. Bcthc, Phys. Rev. Lett. 56, 1305 (1986). 

24. W. C. Haxton, Phys. Rev. Lett. 57, 1271 (1986). 

25. S. J. Parke, Phys. Rev. Lett. 57, 1275 (1986). 

26. C. S. Lim and W. J. Marciano, Phys. Rev. D 37, 1368 (1988). 

27. E. Kh. Akhmcdov, Sov. J. Nucl. Phys. 48, 382 (1988). 

28. G. Barr, T.K. Gaisser, and T. Stanev, Phys. Rev. D 39, 3532 (1989); 
T.K. Gaisser, T. Stanev, and G. Barr, ibid. 38, 85 (1988). 

29. H. Lee and SA. Bludman, Phys. Rev. D 37, 122 (1988); E.V. Bugaev 
and V.A. Naumov, Phys. Lett. B 232, 391 (1989); M. Honda, K, Kasa- 
hara, K. Hidaka, and S. Midorikawa, Phys. Lett. B 248, 193 (1990); M. 
Kawasaki and S. Mizuta, Phys. Rev. D 43, 2900 (1991). 

30. S. Hatakeyama et al. (Kamiokandc collaboration), Phys. Rev. Lett. 81, 
2016 (1998). 

31. M. Ambrosio et al. (MACRO collaboration), Phys. Lett. B 434, 451 
(1998). 

32. A. W. P. Poon, hep-cx/0211013/ 

33. Chooz Collaboration, M. Apollonio et al., Phys. Lett. B 416, 415 (1999). 

34. F. Bochm ct al., Phys. Rev. D 64, 112001 (2001). 

35. C. W. Walter, hep-ex/0211059/ 

36. KamLAND Collaboration, hep-ex/0212021/, to appear in Phys. Rev. 
Lett. 

37. C. Weinheimer, hep-ex/0210050/ 



59 



38. A. Mezzacappa et al., Ap. J. 495, 911 (1998); H.-Th. Janka and E. 
Muller, Astron. Astrophys. 306, 167 (1996); A. Burrows, S. Hayes, and 
B. A. Fryxcll, Ap. J. 450, 830 (1995). 

39. M. Th. Keil, G. G. Raffclt, and H.-T. Janka, astro-ph/0208035/ (sub- 
mitted to Ap. J.) 

40. S. E. Woosley and W. C. Haxton, Nature 334, 45 (1988); S. E. Woosley, 
D. H. Hartmann, R. D. Hoffman, and W. C. Haxton, Ap. J. 356, 272 
(1990). 

41. F. X. Timmcs, S. E. Woosley, and T. A. Weaver, Ap. J. Suppl. 98, 617 
(1995). 

42. Y.-Z. Qian, to appear in Prog. Part. Nucl. Phys. 

43. J. J. Cowan and C. Sneden, astro-ph/0212149/ 

44. S. E. Woosley, J. R. Wilson, G. J. Mathews, R. D. Hoffman, and B. S. 
Meyer, Ap. J. 433, 229 (1994). 

45. G. C. McLaughlin, J.M. Fetter, A.B. Balantckin, and G. M. Fuller, Phys. 
Rev. C 59, 2873 (1999). 

46. W. C. Haxton, K. Langanke, Y.-Z. Qian, and P. Vogel, Phys. Rev. Lett. 
78, 2694 (1997) and Phys. Rev. C 55, 1532 (1997). 

47. G.M. Fuller, RW. Mayle, B.S. Meyer, and J.R. Wilson, Ap. J. 389, 
517 (1992); RW. Mayle, in Supernovae, A.G. Petschek, Ed., (Springer- 
Verlag, New York, 1990). 

48. LSND Collaboration, Phys. Rev. Lett. 75, 2650 (1995); W. C. Louis, 
Nucl. Phys. B 91, 198 (2001); also see K. Eitel (KARMEN Collabo- 
ration), Contribution to the International Nuclear Physics Conference, 
Paris, August 1998. 

49. Y.-Z. Qian ct al., Phys. Rev. Lett. 71, 1965 (1993). 

50. B.S. Meyer et al., Ap. J. 399, 656 (1992)); K. Takahashi,J. Witti, and 
H.-Th. Janka, Astron. Astrophys 286, 857 (1994). 

51. B.S. Meyer, Ap. J. 449, L55 (1995). 

52. G.M. Fuller and B.S. Meyer, Ap. J. 453, 792 (1995). 

53. B. S. Meyer, G. C. McLaughlin, and G. M. Fuller, Phys. Rev. C 58, 
3696 (1998). 



60 



